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ABSTRACT 


A  description  of  the  recording  equipment  and  geological  configuration  at 
the  location  of  three  Soviet  geomagnetic  observatories  which  participated  in 
the  IGY  is  followed  by  a  discussion  of  the  seimple  data  processing  procedures 
for  power  spectral  analysis.  Theory  and  applications  of  time-domain  numerical 
filter  operators  are  presented.  A  discussion  of  the  application  of  the  Tukey- 
Blactanan  method  of  power  spectral  analysis  to  single  time  series  is  followed 
by  its  extension  to  cross-power  spectra,  coherency  and  phase  analysis  of  two 
time  series.  Power  spectra  of  the  horizontal  components  of  the  magnetic  field 
and  telluric  field  are  computed  for  ten-hour  periods  recorded  on  1  and  2 
September  1957  at  the  Soviet  Geomagnetic  Observatory  in  Tbilisi.  Spectral 
analyses  of  the  East-West  telluric  field  component  were  carried  out  for  simul¬ 
taneous  recordings  of  several  ho\irs  duration  made  on  20  September  1957  at 
three  Soviet  stations:  Tbilisi,  Lvov,  and  Ashkhabad.  All  analyses  are  based 
on  microfilm  copies  furnished  by  the  IGY  World  Data  Center  A.  Each  power 
spectrvm  shows  a  frequency  dependence  over  a  frequency  range  of  O.OJ  to  7  cpks 
of  the  form  PqY*”  where  n  varies  from  1.1  to  2.5*  The  coherency  between 
orthogonal  components  of  the  telluric  and  magnetic  fields  is  computed  for  the 
data  from  Tbilisi.  Magnetic  power  spectra  from  USSR,  Canada,  Texas,  and 
Massachusetts  are  plotted  on  a  common  graph  to  show  the  frequency  dependence 
over  the  range  O.OJ  to  400  cpks.  The  power  spectra  for  Tbilisi  are  used  to 
compute  apparent  resistivity  in  accordance  with  a  formula  developed  by 
Cagniard.  It  is  found  that  both  components  of  the  telluric  field  yield  a 
value  of  6  ohm  meters  for  the  electrical  resistivity  of  the  earth.  It  is 
suggested  that  this  value  applies  to  depths  from  70  to  700  kilometers.  There 
is  no  evidence  of  horizontal  anisotropy  of  the  earth  *s  resistivity.  Assuming 
that  the  earth  behaves  as  a  linear  filter  with  the  variation  of  a  horizontal 
component  of  the  earth’s  magnetic  field  as  the  random  input  or  driving  force 
and  the  corresponding  component  of  the  telluric  field  as  the  output  or  result¬ 
ing  function,  the  frequency  response  of  the  earth  is  computed  using  the  results 
of  cross -spectral  analysis.  Resultant  computation  shows  the  earth  to  behave 
as  a  uniform  conductor.  A  comparison  is  made  between  frequencies  at  which 
peaks  occur  in  the  power  spectra  of  the  magnetotelluric  field  and  the  observed 
modes  of  free  oscillation  of  the  earth. 
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INTRODUCTION 


This  dissertation  is  an  investigation  of  the  power  spectra  of  the 
magnetotelliiric  field  of  the  earth  and  some  geophysical  interpretations  of 
this  power  distribution.  The  frequency  range  considered  is  0.03  to  10  cycles 
per  kilosecond,  a  frequency  range  not  previously  investigated  by  power  spectral 
analysis. 

Between  any  two  points  on  the  surface  of  the  earth  there  exists  a  time 
dependent  difference  of  potential.  The  earth,  a  conductor,  provides  a  path 
of  CTirrent  flow.  The  earth  current  or  telluric  current  as  it  is  sometimes 
called,  ranges  in  magnitude  from  zero  to  several  millivolts  per  kilometer  and 
is  easily  measured.  Many  stations  called  earth  current  observatories,  have 
been  set  up  to  make  continuous  recordings  of  telluric  currents  to  provide  a 
basis  for  systematic  study.  Historical  accounts  are  plentiful  and  complete 
and  will  not  be  considered  here  [l,  2,  5]**  Ib  is  sufficient  to  say  that 
earth  currents  exist  and  are  measurable. 

Measurements  of  the  magnetic  field  of  the  earth,  near  the  surface  of  the 
earth,  show  this  field  is  time  dependent.  It  has  been  noted  that  the  time 
variations  in  the  earth  cxirrents  and  in  the  magnetic  field  at  the  surface  of 
the  earth  are  at  times  very  similar  in  appearance  at  the  same  station.  Certain 
characteristics  such  as  world  wide  appearance,  relations  to  simspots,  27-day 
recurrence  and  tendency,  simultaneous  occurrence  of  disturbances,  and  so  forth, 
are  common  to  the  two  fields  at  many  stations.  It  is  generally  believed  that 
this  similarity  is  more  than  chance  and  that  the  two  are  in  some  way  related. 

Since  electrical  measuring  and  recording  instrxaments  are  all  frequency 
sensitive,  the  present  investigation  is  limited  by  the  characteristics  of  the 
recordings  at  hand  for  investigation  and  analysis.  The  recording  instruments 
used  at  the  Soviet  observatories  respond  to  dc  and  the  only  frequency  limita¬ 
tion  in  the  response  occurs  at  frequencies  hi^er  than  those  considered  in 
the  present  investigation. 

Numbers  in  brackets  refer  to  references  in  the  bibliography. 
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The  analysis  here  is  different  from  the  traditional  approach  in  that, 
instead  of  concentrating  on  discrete  events  such  as  storms  and  sudden  com¬ 
mencements,  the  analysis  is  limited  to  carefully  selected  records  which  most 
nearly  approximate  stationary  time  series.  This  selection  is  important  since 
the  data  are  subjected  to  spectral  analysis.  Time-domain  numerical  filters 
are  used  to  implement  the  spectral  analysis. 

The  data,  recording  equipment  and  recording  stations  are  discussed  in 
Chapter  I.  The  various  theoretical  aspects  of  spectral  analysis  and  filter¬ 
ing  are  considered  in  Chapters  II  through  V.  The  power  spectra  are  presented 
in  Chapter  VI.  The  geophysical  interpretations  are  found  in  Chapter  VII. 


CHAPTER  I 


THE  DATA  AND  DATA  PROCESSING 


Source  of  Data 

A  planned  effort  by  the  scientists  of  the  USSR  for  recording  of  the 
variations  in  the  earth's  electromagnetic  field  has  resvilted  in  a  series  of 
zonal  and  meridional  chains  of  recording  stations  extending  across  the  USSR 
and  from  the  Arctic  to  the  Antarctic.  A  discussion  of  this  effort  including 
description  of  the  equipment  is  available  as  part  of  the  Soviet  contribution 
to  the  International  Geophysical  Year  [4]. 

The  data  on  vhich  this  dissertation  is  based  were  obtained  from  the  IGY 
World  Data  Center,  Washington,  D.  C.  Several  rolls  of  microfilm  copies  of 
normal-run  and  rapid-rtin  telliarograms  and  magnetograms  selected  from  records 
deposited  at  the  data  center  by  the  USSR  were  obtained  through  the  Library 
of  The  Uhlversity  of  Texas.  Records  were  obtained  for  the  following  stations: 


USSR 

A020 

Cape  Chelyushin 

USSR 

B145 

Lvov 

USSR 

CI26 

Ashkhabad 

USSR 

C364 

Tbilisi 

USSR 

A9T6 

Oasis 

The  records  made  at  Cape  Chelyushin  and  Oasis  were  unsuitable  for  the 
purpose  at  hand.  The  remaining  stations,  Lvov,  Ashkhabad,  and  Tbilisi  are 
three  of  the  six  stations  referred  to  as  the  "third  zonal  chain"  by 
V.  Troitskaya  [5].  They  are  located  within  a  geographic  band  extending  from 
38“N  to  49"N  latitude  (see  Figure  l). 

Description  of  Observatories  and  Equipment 

The  information  presented  here  is  obtained  from  reports  to  the  Committee 
on  Observatories  of  the  International  Association  of  Geomagnetism  and 
Aeroncany  [4]. 

Tbillol 

Tbilisi  Observatory  is  located  near  the  town  of  Dusheti,  57  kilometers 
northwest  of  Tbilisi  in  the  Georgian  Soviet  Republic,  USSR.  The  geographic 
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location  is  45*5*28”N,  44*42*16'’E,  the  geomagnetic  location  is  36®7'N, 
122®1*E,  and  the  altitude  is  981*8  meters  above  sea  level.  The  recordings 
are  marked  with  Greenwich  Mean  Time  (GMT),  which  differs  from  local  time  by 
two  hours,  58  minutes,  49  seconds.  The  observatory  is  located  in  the  eastern 
suburb  of  the  town  on  a  small  shrub-covered  hill.  An  electric  railway 
operates  approximately  26  kilometers  from  the  observatory,  but  as  far  as  can 
be  determined,  it  has  no  measurable  local  effect.  It  is  believed  that  no 
artificial  disturbances  are  present  which  affect  the  recordings. 

The  following  quotation  [4]  gives  the  geological  configuration  of  the 
recording  area: 

The  northern  rim  of  the  Gorl-Mukhransk  depression,  on  which  the 
Dusheti  Geophysical  Observatory  is  located,  Is  composed  of  Tertiary 
and  Quaternary  deposits  of  great  thickness:  sandstones,  carbonate 
block  conglomerates  of  this  upper  Cretaceous,  clayey  shales,  clays, 
solidly  cemented  conglomerates  of  the  middle  Jurassic,  etc.  The 
Tertiary  rocks  are  covered  with  a  mantle  of  valley  alluvixam  up  to 
150  m  in  thickness.  Electric  prospecting  carried  out  in  the  environs 
of  the  station  shows  that  the  surface  deposits  have  specific  resist¬ 
ances  that  do  not  exceed  50  ohms/ra;  low  values  of  the  specific 
resistance  up  to  10  ohras/m  are  characteristic  of  clayey  deposits  of 
the  northwestern  part  of  the  station  site.  Studies  conducted  in  the 
area  of  the  station  by  the  method  of  telluric  currents  and  vertical 
electric  soundings  have  established  the  absence  of  anomalous  fields 
near  the  station  and  a  sufficiently  homogeneous  geoelectric  structure 
of  low  resistance  (conglomerates  and  sandstones  of  the  Dusheti  forma¬ 
tion).  The  geomagnetic  field  of  the  environs  of  the  Dusheti 
Observatory  is  well  studied.  Magnetic  anomalies  are  absent  over  a 
considerable  area  adjoining  the  station. 

The  magnetographs  at  Tbilisi  consist  of  two  sets  of  the  raediiam  model 
Eshenhagen  [6]  type  (imported  from  Germany  in  1930 ).  One  of  these  is  used  as 
the  recorder  and  the  other  for  control  purposes.  Absolute  measurements  of  D, 
H,  and  Z  are  made  on  the  average  of  six  times  per  month  at  the  observatory. 
!Bie  time  marks  for  the  recording  set  are  made  by  use  of  a  Gaslu  contact  clock 
with  an  average  accuracy  of  ±5  seconds - 

The  variations  in  the  telluric  field  are  recorded  by  automatic  photo- 
galvanographs.  The  mirror  galvanometers  installed  in  the  photo-galvanographs 
have  a  current  sensitivity  of  the  order  lO”  to  lO'^  amperes.  Lead  elec¬ 
trodes  are  groiinded  in  clayey  deposits  (at  a  depth  of  two  meters)  at  the  ends 
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of  two  mutually  perpendicular  one  kilometer  lines.  The  lines  are  oriented 
latitudinally  and  longitudinally  to  measure  the  East -West  and  North-South 
components  of  the  telluric  fields,  respectively. 

Ashkhabad 

Ashkhabad  Observatory  is  located  near  the  village  of  Vannovsky,  21 
kilometers  west  of  Ashkhabad  in  the  Turhman  Soviet  Republic,  USSR.  The  geo¬ 
graphic  location  is  3T**57*0*N,  58®06.5’E,  the  geomagnetic  location  is 
155. 1®E,  and  the  altitude  is  570  meters  above  sea  level.  The  recordings  are 
marked  with  Greenwich  Mean  Time  (GMT)  which  differs  from  local  time  by 
3  hours,  52  minutes,  20  seconds.  The  following  quotation  [4]  gives  the  geo¬ 
logical  configuration  on  the  recording  area: 

The  station  is  located  on  the  front  range  of  the  Kopet-Dag,  in 
the  Firyuzinka  River  valley.  The  valley  is  composed  of  a  sedi¬ 
mentary  Cretaceous  deposits  (sandstones,  argillites,  aleurilites) , 
and  is  covered  with  several  meters  of  proluvium.  Tlie  rocks  have 
a  preferential  SSE  dip  and  NNW  strike.  The  magnetic  pavilions  are 
situated  at  the  northern  foot  of  the  mountain  100  m  from  the  geo¬ 
physical  station.  The  earth -ciirrent  station  is  about  5OO  meters 
above  sea  level  with  maximum  relative  rises  of  100-200  meters/km. 

There  are  no  artificial  disturbances  near  the  observatory. 

Barth  current  recording  equipment  is  on  hand  for  22inm/hour,  90  mm/hour, 
and  30  mm/min.  The  22  ram/hour  recording  is  used  here.  Short  lines  in  the 
form  of  a  cross  of  370  m  length  connect  lead  ribbon  electrodes  (13  in  long  and 
3»5  ni  wide)  laid  at  a  depth  of  2  to  2.5  meters  near  the  outcroppings  of  base¬ 
ment  rocks.  The  lines  are  placed  5O-60  cm  underground.  The  electrodes  are 
spread  out  in  a  zigzag  fashion.  The  interelectrode  resistances  reach 
150  ohms  on  the  North-South  line  and  130  ohms  on  the  East -West  line. 

Lvov 

Lvov  station  is  located  in  the  Trans -Carpathian  region  of  the  Ukraine 
SSR,  10  km  southeast  of  the  town  of  Uzhgorod.  The  geographic  location  is 
48®53*2*N,  22°25.8’E,  the  geomagnetic  location  is  47.0*N,  104.2®E,  and  the 
altitude  is  I50  meters  above  sea  level. 

The  following  description  of  the  geological  configuration  is  quoted  from 
the  International  Association  of  Geomagnetism  and  Aeronomy  [4]: 


-  6  - 


The  locality  is  mountainous  with  relative  eminences  up  to  100  m. 

The  region  is  composed,  basically,  of  volcanogenic  and  sedimentary 
rocks  of  neogen  age.  Under  an  upper  layer  of  sandy  loam  and  loess 
of  thickness  5-10  m  and  specific  resistance  5-5  ohm/ra,  is  a  volcano¬ 
genic  thickness  of  effusive  rocks  with  underlying  tuffs,  tuff  breccia 
and  tuffites  with  an  aggregate  thickness  of  about  500  meters  and  a 
8X)ecific  resistance  up  to  200  ohm/m.  The  volcanogenic  layer  is 
distinguished  by  pronounced  anisotropy.  The  lower  lying  strata  of 
coal-bearing  lagoon-continental  deposits  of  thickness  about  700  m 
has  an  elevated  electrical  conductivity.  At  a  depth  of  15OO  m  is  a 
layer  of  high  resistance  that  is  used  as  the  reference  standard  in 
electrical  mapping. 

Equipment  is  available  for  90  mm/hour,  30  mm/min,  and  22  mm/hour  earth 
current  recordings.  The  90  mm/hour  records  are  used  here. 

Selection  of  Data  for  Processing 

Microfilm  copies  of  the  data  were  obtained  for  September,  1957,  but  no 
attempt  was  made  to  process  all  of  the  data.  The  time  and  expense  of  power 
spectral  analysis  makes  it  desirable  to  plan  the  work  carefully  to  insiire  that 
the  goals  of  the  study  can  be  achieved  with  a  minimiam  of  data  processing.  In 
addition  to  these  considerations  are  the  questions  of  q\xality  in  the  records. 
Records  with  sudden  commencements,  storms,  frequent  off -paper  sections,  and 
very  quiet  records  were  rejected.  Care  was  taken  to  select  records  which  were 
moderately  active  but  appeared  stationary.  Figures  2  to  7  are  reproductions 
of  typical  records  used  in  the  study. 

Pre-Machine  Processing  of  Data 

The  most  suitable  records  were  photographically  enlarged  from  the  35  rom 
microfilm  to  approximately  17  x  22  inches.  Portions  of  the  records  near  the 
edges  of  the  enlargement  were  not  used  since  lens  distortion  during  the 
enlargement  process  was  evident.  The  records  were  then  digitized  by  hand. 

The  accxiracy  of  digitized  data  is  H  gamma  on  the  magnetic  field  records,  ±1 
minute  on  the  declination  records,  and  ±0.01  mv/km  on  the  telluric  records. 

The  accuracy  here  is  considerably  greater  than  necessary  since  the  statisti¬ 
cal  accuracy  of  the  computed  spectral  estimates  is  the  controlling  factor  [7]* 
The  selection  of  the  sampling  interval  was  based  on  the  appearance  of 
the  record.  Approximate  spectral  estimates  were  made  by  visml  inspection 
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FIGURE  3 

MAGNETOGRAM  FOR  TBILISI  2-3  SEPTEMBER  1957 
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of  the  record  to  assure  selection  of  a  sampling  interval  fine  enough  to 
prevent  distortion  of  the  spectrum  throu^  aliasing  (see  Chapter  II). 

The  magnetograms  contained  records  of  the  horizontal  component  of  the 
magnetic  field  (H),  the  vertical  component  of  the  magnetic  field  (Z),  and 
the  declination  (d).  The  H  and  D  records  were  used  to  compute  the  North- 
South  and  East -West  components  of  the  horizontal  magnetic  field  through  the 
relations 

N-S  Magnetic  =  H  cos  D,  and 

E-W  Magnetic  =  H  sin  D. 

Machine  Processing 

Figure  8  is  a  flow  chart  showing  the  steps  in  the  machine  computation  of 
the  power  spectrum  of  a  single  time  series.  In  cases  where  simultaneous  mag¬ 
netic  and  telluric  data  are  processed^  the  program  is  considerably  more 
complicated  since  various  combinations  of  cross-correlation,  cross-power, 
coherency  and  phase  angles  along  with  apparent  resistivity,  apparent  conduc¬ 
tivity,  and  the  frequency  response  of  the  earth  are  computed. 

In  the  early  stages  of  the  work,  computation  was  performed  on  the  IBM  65O 
at  the  Computation  Center  of  The  Iftiiversity  of  Texas.  Subsequently,  the 
IBM  6^  was  replaced  by  a  CDC  160U  machine.  The  bulk  of  the  ccoputation  was 
perfonied  on  the  new  machine. 


FIGURE  8 

MACHINE  COMPUTATION  FLOW  CHART  FOR 
A  SINGLE  TIME  SERIES 
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CHAPTER  II 


SAMPLE  DATA  CONSIDERATIONS 

A  continuous  function  of  time  is  completely  determined  by  its  values  at 
equally  spaced  intervals  provided  that  the  continuous  function  contains  no 
frequencies  higher  than,  say,  W  cycles  per  second,  and  the  ordinates  are  given 
at  points  spaced  l/2¥  seconds  apart,  the  series  of  values  extending  for  all 
time  [8].  Here  the  analysis  is  to  be  based  on  sampled  values  obtained  from 
continuous  records  that  are  not  infinite  in  extent  and  are  not  band  limited. 
Analysis  based  on  finite  amounts  of  data  is  common  to  statistical  work  and 
will  be  considered  in  due  course.  The  fact  that  the  original  functions  are 
not  band  limited  will  be  discussed  in  connection  with  the  problem  of  aliasing. 

Aliasing 

An  important  consideration  in  the  computation  of  power  spectra  from 
equi -spaced  sample  data  is  the  problem  of  aliasing.  Figure  9  shows  portions 
of  two  sine  waves  of  different  frequencies.  Attention  is  focused  on  the 
equally  spaced  sample  data  values  which  are  the  same  for  each  sine  wave  even 
thou£^  the  waves  differ  in  frequency.  Thus,  given  only  the  sample  values,  a 
sine  wave  of  a  given  frequency  may  be  confused  with  a  sine  wave  of  hippier 
frequency.  Specifically,  if  a  harmonic  time  function  x(t)  is  sampled  at 
equally  spaced  time  Intervals  At,  then  a  frequency 

f 

N  2At  ’ 

called  the  Nyquist  or  folding  frequency  [9],  exists  such  that  the  functions 
with  frequencies 

f  +  n  f N  ,  for  n  =  0,2,4,-**  , 

are  not  distinguishable.  For  example,  consider  the  function 

x(f,t)  =  sin  2itft  , 

sampled  at  kAt  (k  =  0,  1,  2,**«)  so  that  the  sampled  values  of  x  are 
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x(f,l^t)  =  sin  2itfk 
■ov  consider  the  function 

x(f  ±  n  *  sin  2*(f  +  n  fjj)t 


'k  ■  St  > 
n  .  0,  2,  4,  6,"-, 

ssaqpled  at  equl -spaced  Intervals  At  so  that  the  sampled  values  are 
x(f  ±  n  fjj,  l^t)  =  sin  2K(f  ±  n  fjj)kAt 

Use  of  simple  trigonometric  Identities  reveals  that 
x(f,  kAt)  =  x(f  ±  n  fjj,  kAt) 

Obviously,  then,  power  contributed  to  a  power  spectrum  at  a  given  frequency  f 
cannot  be  distinguished  from  powers  contributed  by  frequencies  f  ±  n  f^^. 

This  translation  of  frequencies  is  known  as  aliasing  [7l*  If  the  data 
actually  contain  power  at  frequencies  greater  than  f^^,  this  power  will  be 
"folded  back"  into  the  principal  band  which  extends  from  0  to  f^^.  Power  that 
is  folded  back  results  in  a  distortion  of  the  true  power  spectrum  in  the 
principal  band. 

To  make  the  effect  of  aliasing  negligible  it  is  necesscury  to  select  a 
sampling  interval  "small  enough"  to  place  the  Nyqulst  frequency  beyond  all 
significant  power  contributions. 

Frequency  Resolution 

Questions  of  frequency  resolution  have  to  do  with  the  number  of  power 
spectral  estimates  per  unit  frequency  and  to  what  extent  a  particular  estimate 
has  been  affected  by  power  at  other  frequencies. 

Leaving  aside  for  the  moment  the  possibility  of  filtering,  the  selection 
of  the  length  of  the  sampling  interval.  At,  is  determined  by  the  relative 
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amount  of  power  occurring  at  frequencies  higher  than  the  Nyquist  frequency 
associated  with  At.  The  sampling  Interval  must  he  selected  small  enough  to 
bring  the  aliased  power,  if  any,  to  a  tolerable  minimian.  That  is,  the  con¬ 
tribution  due  to  aliased  power  should  be  less  than  the  statistical  accviracy 
of  an  estimate  in  the  principal  band.  Once  the  length  of  sampling  interval 
has  been  established  and  the  autocorrelation  is  computed  for  m  lags  (see 
equation  31)^  the  power  spectrum  is  computed  at  m  +  1  equally  spaced  fre¬ 
quency  values  in  a  bemd  extending  from  zero  to  the  Nyquist  frequency  [9]. 

For  a  fixed  record  length  a  compranlse  must  be  reached  between  frequency 
resolution  and  the  stability  of  power  spectral  estimation.  Attempts  to 
Increase  the  frequency  resolution  by  increasing  the  number  of  lags  are  met 
by  a  decrease  in  the  stability,  since  the  stability  is  a  function  of  the 
number  of  lags  [7].  However,  an  increase  in  the  frequency  resolution  can  be 
effected  without  a  corresponding  decrease  in  stability  by  increasing  the 
record  length.  This  procedure  may  lead  to  a  prohibitive  amoimt  of  additional 
labor. 

The  refined  power  spectral  estimate  at  a  particular  frequency  is 
computed  using  a  weighted  average  of  the  raw  estimate  at  that  frequency  and 
the  raw  estimates  immediately  adjacent  to  that  particular  frequency  (see 
equations  3^)  35,  and  36). 

Purpose  of  Filtering 

In  the  spectral  computations  considered  here,  a  high-pass  filter  is  used 
to  prewhiten  (specifically,  flatten)  the  spectrum,  and  a  low-pass  filter  is 
used  to  reshape  the  spectrum  prior  to  decimation  (see  next  page). 

A  common  feature  of  the  power  spectra  computed  and  plotted  here  is  the 
rapid  decrease  of  power  with  increasing  frequency.  Tukey  points  out  how  the 
computation  of  estimates  on  a  spectrum  with  large  slope  leads  to  errors.  The 
purpose  of  high-pass  filtering  is  to  make  the  spectrum  flat  before  the  power 
spectral  computation.  The  filter  is  removed  after  computation  to  restore 
the  spectrum  to  its  true  shape. 

When  it  is  of  interest  to  examine  the  low  frequency  portion  of  the 
Hljrquist  band  in  greater  detail,  that  is,  compute  spectral  estimates  at  closer 
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frequency  intervals,  two  ways  of  accomplishing  this  are  available.  As  has 
already  been  mentioned,  the  number  of  lags  could  be  increased.  This  requires 
a  large  amount  of  additional  labor  and  for  a  fixed  length  leads  to  increased 
instability  in  the  estimates.  The  other  approach  is  to  apply  a  low -pass  fil¬ 
ter  to  the  original  data,  with  sampling  interval  At,  and  then  compute  the 
power  spectrum  with  respect  to  a  new  Nyquist  frequency  f^.  The  process  of 
converting  the  sampling  interval  to  a  new  sampling  interval  which  is  an  inte¬ 
ger  multiple  of  the  original  interval  is  known  as  decimation. 

If  it  is  possible  to  compute  an  undistorted  spectrum  after  decimation, 
then  the  original  At  was  smaller  than  necessary.  In  general,  decimation  must 
be  preceded  by  low-pass  filtering.  Figure  10  shows  an  example  of  a  power 
spectrum  before  and  after  low-pass  filtering. 

Confidence  Limits 

The  values  of  power  spectral  density  are  computed  from  finite  amounts 
of  data  and  therefore  are  statistical  estimates  of  the  true  values.  The 
Tukey  method  [7 ]  provides  confidence  intervals  by  which  the  variability  of 
the  estimates  may  be  computed.  Confidence  intervals  are  a  familiar  feature 
of  statistical  analysis  [lO]  and  will  not  be  discussed  here.  However,  it  is 
important  to  note  that  the  basis  for  the  Tukey  values  is  the  assumption  that 
the  process  is  Gaussian  and  the  distribution  of  the  spectral  estimates 
follows  the  so-called  "chi-square"  distribution. 

All  of  the  power  spectra  computed  here  have  eighty  percent  confidence 
limits  indicated  on  the  graphs. 
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POWER  SPECTRUM  BEFORE  LjOW-MSS  FILTERING 


POWER  SPECTRUM  AFTER  LOW-FASS  nUERiNG 


FIGURE  10 

EXAMPLE  OF  A  POWER  SPECTRUM  BEFORE  AND  AFTER 
LOW-PASS  FILTERING 
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CHAPTER  III 


TIME-DOMAIN  FILTERS 

Figure  11  shows  a  block  diagram  of  a  stationary  linear  passive  electric 
wave  filter.  A  four  terminal  electrical  network  is  said  to  be  a  stationary 
linear  passive  filter  if  (l)  the  properties  of  the  network  do  not  change  with 
time,  (2)  the  output  voltage  remains  at  zero  xmtil  a  zero  input  voltage  changes 
to  a  non-zero  value,  (3)  a  linear  combination  of  input  voltages  results  in  an 
output  voltage  that  is  the  sum  of  the  effects  due  to  each  input  voltage 
8ex>arately,  and  (4)  there  are  no  energy  sources  within  the  circuit.  For  a  sta¬ 
tionary  linear  passive  filter,  the  output  voltage,  expressed  as  a  function 
of  time,  may  be  obtained  from  the  input  through  some  relating  operator. 

Here  the  interest  is  in  linear  operators  which  mathematically  perform  the 
same  operation  as  an  electrical  filter,  but  are  more  general  in  that  they 
are  not  limited  by  the  realizability  conditions  imposed  on  physical  filters. 
Both  the  impulse  response  and  the  frequency  response  are  characteristic  of 
the  filter  for  particular  inputs.  That  is,  the  impulse  response  is  the 
output  for  a  unit  impulse  (Dirac  delta  at  t  =  O)  input,  whereas  the  fre¬ 
quency  response  is  the  relation  between  the  output  and  input  for  a  sinusoidal 
input.  In  particular  the  impulse  response  relates  the  output,  H(t),  to  the 
input,  G(t),  through  the  relation  [ill 

P  00 

H(t)  =  /  G(t  -  t)  W(T)dT  ,  (1) 

where  G(t)  is  the  input  at  time  t 

H(t)  is  the  output  at  time  t,  and 
W(t)  is  the  impulse  response. 

TOie  0  extends  the  lower  limit  to  negative  values  since  the  Dirac  delta  input 
at  t  =  0  extends  "across  zero."  The  impulse  response  shows  the  "memory” 
proi)erty  of  the  filter  in  that  it  determines  how  the  past  of  the  input 
effects  the  output  at  time  t. 

If  the  input  is  a  unit  amplitude  sinusoidal  function,  say. 
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FIGURE  li 

FOUR  TERMINAL  ELECTRICAL  FILTER 


-  2L  - 


(8) 


0(t)  = 

and  then  one  uses  equation  (l)  to  find  the  output,  one  has 

H(t)  -  r  W(t)  ,  (3) 

'jQ 

or 

H(t)  -  r  W(t)  (4) 

^O" 

Letting 

Y(f)  =  r  W(t)  ,  (5) 

one  has 

H(t)  =  Y(f)  G(t)  .  (6) 

Y(f),  as  defined  In  eqioatlon  (5)  Is  what  Is  usually  referred  to  as  the  fre¬ 
quency  response  of  the  filter. 

In  this  paper  filtering  will  be  accomplished  by  a  linear  operation 
referred  to  as  time-domain  filtering  [12].  The  filter  Input  consists  of  a 
set  of  numbers  which  represent  the  amplitude  of  the  record  to  be  emalyzed  at 
equally  spaced  points  separated  by  the  time  Interval  At.  The  filtering  opera¬ 
tion  is  carried  out  by  using  the  relation 

K 

v(t)  *  ^  -  kAt)  ,  (7) 

k=0 

where  v(t)  is  the  filter  output  at  t, 

u(t)  is  the  filter  input  at  t,  and 

aj^  are  a  set  of  K  +  1  numerical  coefficients  which  will  be  referred 
to  as  the  filter  coefficients  of  the  linear  operator. 

Equation  (7)  may  be  considered  as  an  approximation  to  equation  (l)  [ij]. 

This  is  acconpllshed  by  replacing  the  Infinite  limit  by  the  finite  number  K, 
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applying  Simpson's  rule  to  the  integrand  (using  At  as  the  spacing  interval), 
and  choosing  the  values  of  W  to  be 


W(lsAt)  =  aj^ 


(8) 


where  k  =  0,  1,  2, • • • ,  K. 

In  his  paper  [13],  Dr.  Horton  discusses  the  mean  square  error  in  this  approxi¬ 
mation  and  a  method  for  calculating  the  numerical  coefficients  to  be  used  to 
obtain  a  desired  frequency  response.  The  fact  that  this  is  an  approximation 
will  be  suppressed  here  since  the  exact  values  of  the  frequency  response  for 
a  corresponding  physical  filter  need  not  be  considered.  The  viewpoint  is 
taken  here  that,  given  a  function  u(t)  defined  at  discrete  equally  spaced 
values  and  given  a  set  of  K  +  1  numbers,  aj^,  a  new  function  v(t)  may  be 
computed  according  to  the  relation  (7). 

If  the  input  function  happens  to  be  a  unit  sinusoidal  function  exp(j2icft), 
then  the  computed  (filtered  or  output)  function  becomes 


K 

v(t)  «  V  e 


J2jtf(t  -  l^t) 


(9) 


or 


V,.)  .  ^  f 

IlD-O 


Letting 


a^e’J2nfl«^M 


(10) 


Y(f) 


K 

I 

k^O 


a,_  e 


■J2xfkAt 


(11) 
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one  has 


v(t)  =  Y(f)  u(t)  .  (12) 

Y(f),  as  expressed  by  equation  (11),  vlU  he  referred  to  as  the  frequency 
response  function.  Since  Y(f)  Is  periodic  of  period  l/dt,  a  selection  of 

1  <  f  <  j.  1 

■  2At  "  2£^t 

for  the  fundamental  range  of  definition  relates  the  cut-off  frequency  of  the 
filter  to  the  Nyquist  frequency  for  a  sampling  interval  At. 

The  two  filters  used  in  this  paper  are  em  l8-term  hl^-pass  filter  [12] 
and  a  9-tem  low -pass  filter  [l?].  The  numerical  coefficients  for  these 
filters  are  given  In  Table  1. 


!nie  frequency  response  Is  CGisputed  (at  N  1  points)  from  the  filter 
coefficients  by  letting 

N 

DfkO 

-  A(f)  -J  B(f)  , 

SO  that 


and 


N 

A(f)  “  ^  *^n  ®°*(2icnfAt)  , 

N 

B(f)  “  ^ 

1^ 


(15) 


(16) 
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TABLE  I 


LINEAR  OPERATORS  USED  TO  FILTER  RAW  SAMPLED  DATA 
Case  A:  High-pass  Filter  (l8  coefficients) 
Case  B:  Low -pass  Filter  (9  coefficients) 


1 

Case  A 

Case  B 

0 

49.2 

0.0078 

1 

-100.0 

0.0508 

2 

51.5 

0.1266 

5 

72.8 

0.1977 

4 

-13.4 

0.2302 

5 

-33.2 

0.1977 

6 

-17.9 

0.1266 

7 

3.3 

0.0508 

8 

9-9 

0.0078 

9 

3.8 

10 

-2.7 

11 

-3.8 

12 

-1.0 

15 

1.2 

14 

1.2 

15 

0.1 

16 

-0.6 

17 

-0.4 
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Recalling  that  the  fundamental  period  Is  defined  for  the  frequency  range 

0  $  |f|  ■  l/2At  It  Is  convenient  to  set 

f 

N  2^t  ' 

(17) 

and 

f  =  ^  f  „ 
q  m  N 

(18) 

where  fj^  Is  the  Nyqulst  frequency.  Equation  (l8) 
tlone  (15)  and  (l6)  to  give 

Is  substituted  Into  equa 

> 

II 

\  f  5^)  > 

(19) 

and 

N 

B(f,)  =  y 

(20) 

iW5 


where  q  =  0,  1,  2,  5>  “  ’ >  m* 

The  frequency  response  is  then  computed  from 


|Y(f)l  =  1  Aitf  +  B(f)^  1 

(21) 

and  the  phase  delay  Is  given  by 

9(f)  *  arctan 

(22) 
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Figure  12  shows  a  plot  of  the  filter  coefficients  versus  time  for  the 
high-pass  filter.  Case  A  of  Table  1.  Figure  IJ  shows  a  plot  of  the  filter 
coefficients  versus  time  for  the  low-pass  filter,  Case  B  of  Table  1. 

Figure  l4  shows  the  frequency  response  of  the  hlg^-pass  and  low-pass 
filters . 

Figure  1^  shows  the  phase  delay  for  the  high-pass  filter.  The  phase 
delay  for  the  low-pass  filter  is  a  stralj^t  line  since  the  filter  is 
syoMtrlc  and  is  not  plotted  here. 
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FILTER  COEFRCIENT 


<n/2 

FREQUENCY  IN  CYCLES  PER  KILOSECOND 


FIGURE  14 

THE  FREQUENCY  RESPONSE  OF  LOW-PASS 
AND  HIGH-PASS  FILTERS,  TABLE  I. 
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EDGE  OF  THE  FUNDAMENTAL 
PERIOD  OF  THE  OPERATOR 


DELAY  IN  DEGREES 


4360 


4170 


4180 


440 


Oh 


-40h 


-180 


O  O  ^  D  O 


J _ L 


|fN 

FREQUENCY 


FIGURE  15 

THE  PHASE  DELAY  FOR  HI6H-R(\SS  FILTER, 
CASE  A  OF  TABLE  I 
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CHAPTIR  IV 


STATISTICAL  ANALYSIS  FOR  A  SINGLE  TIME  SIBIES 

The  results  presented  In  this  dissertation  are  obtained  by  the  applica¬ 
tion  of  the  technique  of  power  spectral  analysis  to  the  study  of  the  low 
frequency  magnetotelliirlc  field  of  the  earth.  Computational  procedures  are 
used  that  yield  the  statistical  variance  spectrvon  of  a  time  series.  Other 
names  for  the  resultant  computation  are  power  density  spectrum,  second-degree 
spectrum  and  quadratic  spectrum.  All  of  these  names  refer  to  the  distribution 
of  variance  as  a  function  of  frequency.  It  is  beyond  the  scope  of  this  paper 
to  present  a  formal  derivation  of  the  formulas  used  In  power  spectral  computa¬ 
tion;  many  complete  discussions  of  this  type  are  available  [7,  1^,  l6,  1?]. 
However,  an  outline  of  the  results  will  be  presented  to  acquaint  the  reader 
with  the  basis  for  the  formulas  used.  The  excuse  for  the  use  of  the  engineer¬ 
ing  term  "power"  In  connection  with  statistical  variance  Is  that  the  variance 
for  a  time  series  with  zero  mean  Is  equal  to  the  ac  power  dissipated  In  a  one 
ohm  p\ire  resistance  [9]. 

For  deterministic  functions  such  as  periodic  and  aperiodic  functions  a 
harmonic  analysis  is  usually  carried  out  by  Foiurler  series  analysis  and  by 
Fourier  Integral  analysis,  respectively.  The  discrete  line  spectrum  for  a 
periodic  function  and  the  continuous  spectrum  for  the  aperiodic  function  may 
be  determined  analytically  because  these  deterministic  functions  are  by 
definition  "known  for  all  values  of  time." 

Random  series  [17]  are  a  class  of  functions  which  are  not  deterministic 
and  do  not  lend  themselves  to  the  harmonic  analysis  techniques  discussed 
above.  In  this  paper  attention  will  be  restricted  to  those  random  time 
series  which  very  closely  approximate  a  stationary  random  ergodlc  process. 
Functions  which  represent  stationary  random  ergodlc  processes  are  characterized 
by  probability  distributions  that  do  not  change  with  time  and  the  time  average 
over  a  sample  function  of  an  ensemble  of  functions  is  equivalent  to  an  average 
over  the  ensemble  of  functions.  These  functions  are  referred  to  as  stationary 
tine  series.  Records  of  the  earth's  magnetic  said  electric  fields  sore  selected 
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vhlch  are  asevoned  to  approximate  closely  the  above  conditions.  It  should  be 
noted  that  proof  of  the  ergodlc  hypothesis  In  a  practical  case  appears  to  be 
Impossible  [l^]. 

Wiener  [l6]  has  shown  that  the  autocorrelation  function, 
defined  by 

1 

q>^(T)  -  llm  i  /  x(t)  x(t  +  t)  dt  ,  (25) 

**  P400 

for  -oe  <  X  <  +  *, 

where  x(t)  Is  a  stationary  time  series,  has  a  Fourier  transform, 
given  by 

=  f  dx  ,  (24) 

^  ••00 


which  turns  out  to  be  the  power  density  spectrum  of  x(t).  Stated  another  way, 
the  Wiener  theorem  [19]  shows  that  the  power  density  spectrum  of  x(t),  defined 


q)  (f)  =  lim  i  f  x(t) 


bauB  a  Fourier  transform 


■  I 


♦  (f)  e 

XX 


which  turns  out  to  be  the  autocorrelation  function  for  x(t;. 

The  power  sx)ectral  analysis  is  based  on  computational  methods  recom¬ 
mended  by  Blackman  and  Tukey  [?]•  The  method  consists  of  estimating  the 
spectinim  of  a  discrete  time  series  by  computing  numerical  approximations  to 
eqmtions  (25)  and  (24).  Since  equation  (25)  is  an  even  function  of  equa¬ 
tion  (24)  may  be  simplified  to 


P  00 

♦^(f)  =2  J  C08  2Hft  dt 
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An  important  featvire  of  the  Tukey  method  la  that  It  provides  a  means  of  esti¬ 
mating  the  power  spectral  density  from  finite  sets  of  sampled  data  with  a 
knowledge  of  the  accuracy  of  the  estimate.  A  brief  outline  of  the  procedure 
will  be  given,  but  no  attempt  will  be  made  to  reproduce  this  work  [7l  here. 

An  experimental  data  record  of  finite  length  Is  sampled  at  efually 
spaced  Intervals  (At)  yielding  a  discrete  time  series 


Xg,  Xj,*”, 


which  will  be  referred  to  as  the  raw  data.  As  has  already  been  shown 
(Chapter  II)  the  selection  of  the  length  of  the  sampling  Interval  for  a 
particular  record  has  Important  Implications.  To  avoid  trouble  from  allasr 
ing  a  minimal  requirement  on  the  choice  of  At  would  be  such  that  no  more 
than  one  relative  maxima  or  minima  occur  In  any  Interval  At. 

The  average  value  for  the  series. 


N 

<  X(t)  >  =  I  ^ 

n^l 


where  X^  =  X(nAt) 

for  n  =  1,  2,*“,  N, 

Is  subtracted  from  each  X  using  the  relation 


X 


n 


-  <  X(t)  > 


(28) 


(29) 


to  remove  the  dc  or  mean  value  fk'on  the  series. 

The  X^  are  then  filtered  (see  equation  7).  Bach  filtered  value  Is 
computed  by 

K 

*n  "  Z  “k  \  -  l^t  ’  ^50) 

b-0 
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where  x  are  the  filtered  valuee. 

n 

are  the  operator  coefflclente, 
i^t  is  the  sampling  Interval, 

X  are  the  raw  data  from  which  the  average  value  has  been 
subtracted,  and, 

k+1  is  the  number  of  filter  coefficients. 

It  should  be  noted  that  this  method  of  filtering  reduced  slightly  the  number 
of  values  available  for  the  autocorrelation  computation.  Fortunately,  suit* 
able  filters  with  relatively  few  coefficients  were  available. 

The  estimates  of  the  unnormalized  autocorrelation  at  lag  t(*  ndt)  are 
eoeqputed  from  the  filtered  values  by 

N-K-n 

'xx'’*  ■  iP^  I  ’‘i’  ■ 

for  n  =  0,  1,  2, * • • ,  m, 

where  C  (t)  is  the  computed  autocorrelation  at  lag  t(»  ndt), 

JUv 

N  Is  the  niMber  of  raw  data  values  > 

K+1  is  the  number  of  filter  coefficients^ 

m  is  the  number  of  lags^  and 

X  are  the  filtered  values. 

P 

C  (t)  are  often  referred  to  as  the  seunple  serial  products,  mean  lagged 
products,  or  autocorrelation  coefficients.  Figures  l6  and  17  are  autocorre¬ 
lations  computed  from  high-pass  filtered  values  and  from  low-pass  filtered 
values.  Comparison  of  Figures  l6  and  1?  shows  the  obvious  effects  of  filter¬ 
ing  and  decimating. 

The  raw  spectral  density  estimates  are  cominited  by  applying  a  discrete 
finite  cosine  series  transform  to  the  mean  lagged  products  according  to  the 
formula 
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I  LA6  •  72  SECONDS 


■flOjOOO' 


FIGURE  16 


THE  AUTOCORR|^AT)QN_FOR 


DATA  FROM 


TBILISI  I  SErrtMBER  1957 


^  .  HI6H-RASS  FILTERED 

WtST  MAGNETIC  FIELD  AT 


AUTOCORRELATION 


FIGURE  17 

THE  AUTOCORRELATION  FOR  LOW- PASS  FILTERED 
DATA  FROM  EAST-WEST  MAONETIC  FCLO  AT 
TMLtSI  I SEPTEMKR 1967 
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— A 

^(r)  -  Cj^(O)  ♦  2  ^  (^)  *  eo*(iir)  M,(3i) 


for  r  =0,  1,  2, •••,  m. 

It  is  necessazy  to  canpute  only  for  r  in  the  above 

is  qiBMtric  in  r.  The  firefMBey  correqopilni  to  r  Is 


ranfls  sines  Q„(r) 


!nie  systematic  statistical  errors  resulting  ftrom  use  of  a  finite  length 
of  data  appear  in  the  values  of  Ihe  Tukey  technique  to  obtain  improved 

spectral  estimates  involves  a  smoothing  or  refining  operation  performed  on 
the  raw  estimates.  Ibe  specific  smoothing  operation  used  here  is  referred 
to  as  banning  [7].  l!he  refined  spectral  estimates,  are  coaputed  using 
the  relations 

P^(0)  =  0.5  [Q^(0)  +  »  (3‘^) 

P^(r)  -  0.5  ^(r)  +  0.25  (55) 


P^(m)  =  0.5  +  ^(a-l)l#  (36) 

where  r  =  1,  m-1. 

Finally  the  filter  is  removed  by  dividing  each  refined  power  spectral 
estimate  by  the  square  of  the  frequency  response  of  the  filter  at  that  fre¬ 
quency  (see  equation  21).  That  is, 

P'  (f) 

P„(f)  -  ,g  (37) 

lAere  P^(f)  is  the  refined  spectral  estimate  at  f  after  filter  removal, 
P^(f)  is  the  refined  spectral  estimate  at  f,  and 


Y(f)  Is  the  filter  response  at  f. 

It  Is  iiiportant  to  note  that  the  spectral  estimate  as  expressed  In 
equation  (37)  Is  a  function  of  ftrefuaney  only  and  no  InfonMtloa  tStanb  fkaaa 
la  ooBtalaad  tharaln. 
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CHAFTBl  V 


STATI9IICAL  ANALYSIS  FOR  TWO  TIME  SKIES 

Tte  ■tatlstlcal  MAlysls  for  a  sincla  tine  aerlaa  is  esctondad  to  tba 
at  tvo  time  series  by  defining  the  cross -correlation  ftCMtlona  (l8] 

1 

lln  =  /  x(t)  y(t  ♦t)  4t  ,  (36) 

9*^^  J-T/2 


nT/2 

liB  i  /  y(t)  x(t  ♦  t)  4t  ,  (39) 

^  ''-T/2 

time  series  X(t)  and  T(t).  The  correapopdd.ng  eroas- 
are  given  by 

f  dx  ,  (llO) 

f  dr  ,  (41) 

which  are  the  Fourier  transforms  of  equations  (jB)  and  (39) >  respectively. 

Again  the  method  consists  of  computing  estimates  of  equation  (4o)  which  are 
numerical  approximations. 

TWO  experimental  records  of  finite  length  are  sampled  at  equally  spaced 
intervals  (dt)  yielding  the  discrete  time  series 

X^,  Xg,  Xj,***,  Xjj  ,  and 

Yg,  Y3,-“,Yn  . 

The  average  value  is  removed  from  each  series  and  then  each  series  is  subjected 
to  the  filtering  operatlm  (see  equations  29  and  y)), 
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■ 


where 


<  X  <  +  “ 


for  the  stationary 
spectral  densities 


♦  (f)  = 
yx'  ' 


nie  cross-correlation  functions,  equations  (33)  and  (39) >  Are  not 
•ysasetrlc  functions  of  t  and  changing  th*  order  of  x  and  y  changes  the  value. 
iMMvar,  and  cp^  are  related  by 

Furthermore,  unlike  the  autocorrelation  functions,  the  cross -correlation 
functions  are  not  "even"  functions  of  t.  This  means  that  the  Fourier  transfoza 
vlll  have  both  real  and  Imaginary  parts,  nie  estimates  of  the  unnormalized 
eroas-eorrelatlon  at  lag  ±t(t  »  ±n6t)  are  ccsqputed  from  the  filter  values  by 

N-K-n 

P-1 

and 

N-K-n 

prt, 

for  n  =  0,  1,  2, * • • ,  m, 

vhere  C  ( t)  is  the  positive  shift  computed  autocorrelation  at  lag  t, 

C^(-t)  is  the  negative  shift  computed  autocorrelation  at  lag  t, 

N  is  the  number  of  raw  data  values, 

K+1  Is  the  number  of  filter  coefficients, 

m  is  the  number  of  lags,  and,  euid  are  the  filtered  values. 

By  separating  the  odd  and  even  contributions  to  the  cross -correlations, 
the  Fourier  transformation  can  be  carried  out  by  t-niring  the  sina  and  cosine 
transformations  separately.  That  Is,  defining 

I^(T)  ■i[o^(T)  .C^(-t)]  ,  (45) 
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eroM'CpMtral  density  is  computed  by 

^(r)  -  »•  <J^(r)  +  J  Q^(r)|  , 

in  vhleh 


[vH-  1  (g) 


+  R^*“  («)  cos(rit)  j  At  , 


r 

2 

y 

ROdd 

Li 

xy 

J 

-odd 

n 

xy 

for  r  =0,  1,  m, 

vfaere  the  firequency  corresponding  to  r  is 
^  r 


The  real  and  imaginary  parts  of  are  each  refined  by  banning  (see 

equations  and  ^6)  and  the  filter  is  removed  from  each  (see  equation  37)* 

The  resultant  cross-spectral  density  is  written  as 


P  (f)  =  Re[p  (f)l  +  J  I  [p  (f) 
xy'  ^  xy'  '  J  *'  m  I  xy'  * 
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or 

irh«r« 


xy' 


and 


xy 


m  xy' 


)Tl  21 

Jj 


(sit) 


ISie  argiment  of  the  cross -spectral  density  6  (f),  represents  the  phase 

xy 

difference  hetveen  the  power  spectral  density  of  the  X  and  Y  record  at  fre¬ 
quency  f.  Specifically, 

B^{t)  =  e^it)  -  yf)  ,  (55) 

where 


6  (f)  is  the  phase  difference  at  f, 

xy 

©^(f)  is  the  phase  of  the  f  component  of  x  record,  and 
©  (f)  is  the  phase  of  the  f  component  of  y  record. 

y 

Wie  ratio  of  the  magnitude  of  the  cross-power  density  to  the  square  root  of 
the  product  of  the  (auto-)  power  densities  of  the  x  suid  y  records  is  a  mea¬ 
sure  of  the  coherency  between  the  records  [19]»  Hiat  is,  the  coherency 
coefficient  is  defined  by 


Coh^(f) 


(5€) 
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Vrai  tlM  relation  [12] 


it  l8  clear  that  P^(f)  1b  different  from  zero  only  when  common  frequencies 
xy 

exist  in  the  X  and  Y  records.  For  perfectly  coherent  records  the  equality 
prevails  and  the  coherency  is  equal  to  unity.  In  practice^  when  finite  amounts 
of  data  are  used,  the  conputed  values  of  the  coherency  nay  raafe  to  values 
greater  than  one. 
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CHAPTER  VI 


PCWER  SPECTRA  OF  THE  MAGNETOTELLURIC  FIELD 

The  previous  chapters  present  the  methods  of  data  selection,  data 
processing,  and  the  practical  formulas  of  power  spectral  analysis.  This 
chapter  contains  the  computational  results.  Table  2  shows  a  summary  of  the 
data  analyzed  that  are  comimted  and  the  computational  parameters  associated 
with  each  sample.  Geophysical  Interpretations  of  these  computed  results 
are  found  In  Chapter  VII. 

East-West  Telluric  Current  at  Three  Russian  Stations 

Figures  l8,  19,  and  20  are  graphs  of  the  power  sijectra  of  the  East -West 
component  of  the  telliaric  ciorrent  computed  from  several  hoxors  of  rapid-run 
tellurograms  recorded  simultaneously  at  Ashkhabad,  Lvov,  and  Tbilisi  on 
20  September  1957.  The  So^t  confidence  intervals  for  the  decimated  (x)  and 
vindecimated  (o)  spectra  are  indicated  in  each  graph.  The  values  from  which 
the  graphs  are  plotted  are  given  in  Tables  5  through  8.  All  three  graphs 
show  a  characteristic  linear  trend  superposed  by  a  series  of  small  spectral 
peaks.  In  each  graph  a  small  amount  of  aliasing  is  evident  near  the  Nyqulst 
frequency  of  the  spectra  computed  from  low -pass  filtered  data.  If  each  graph 
is  approximated  by  a  straight  line,  it  is  found  that  a  line  of  slope  -1.6 
very  nearly  fits  all  three.  Within  the  limits  of  this  approximation  an 
analytical  relation  of  the  form 

P(f)  =  Pq  f'"  (58) 

where  n  =  -1.6^ 

exists  between  the  power  density  and  frequency  for  the  frequency  range  covered 
by  these  graphs.  The  question  of  an  empirical  formula  suid  the  meaning  of  the 
spectral  peaks  will  be  considered  in  due  course. 

It  is  of  interest  to  compare  the  power  density  at  these  three  stations 
at  a  particular  frequency.  At  1  cpks,  using  Lvov  as  a  reference,  the  power 
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FIGURE  18 

THE  POWER  SPECTRUM  OF  EAST-WEST  TELLURIC 
FIELD  AT  ASHKHABAD  20  SEPTEMBER  1957 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

.141837E+02 

•OOOE+OO 

.000 

.OOOOOOE+00 

2 

.123920E+02 

.278E-04 

600.000 

.866770E+05 

3 

.111016E+02 

.556E-04 

300.000 

.343512E+05 

4 

.102550E+02 

•833E-04 

200.000 

.177925E+05 

5 

.871901E+01 

.lllE-03 

150.000 

. 140715E+05 

6 

.759887E+01 

.139E-03 

120.000 

.990420E+04 

7 

,695464E+01 

.167E-03 

100.000 

.606458E+04 

8 

.653183E+01 

.194E-03 

85.714 

.448033E+04 

9 

.612126E+01 

.222E-03 

75.000 

.337211E+04 

10 

.568163E+01 

.250E-03 

66.667 

.233831E+04 

11 

,535285E+01 

.278E-03 

60.000 

.153280E+04 

12 

•512184E+01 

.306E-03 

54.545 

.898129E+03 

13 

,490455E+01 

.333E-03 

50.000 

.633444E+03 

14 

,490734E+01 

.361E-03 

46.154 

.964231E+03 

15 

.440409E+01 

.389E-03 

42.857 

.168762E+04 

16 

.412987E+01 

.417E-03 

40.000 

.207727E+04 

17 

,399599E+01 

,444E-03 

37.500 

.199364E+04 

18 

.355403E+01 

.472E-03 

35.294 

.211508E+04 

19 

,304744E+01 

•500E-03 

33.333 

.281794E+04 

20 

,293748E+01 

.528E-03 

31.579 

.313296E+04 

21 

.243778E+01 

.556E-03 

30.000 

.218798E+04 

22 

.224384E+01 

.583E-03 

28.571 

. 198563E+n4 

23 

.221846E+01 

•611E-03 

27.273 

.309661E+04 

24 

.192461E+01 

.639E-03 

26.087 

.280845E+04 

25 

,204194E+01 

,667E-03 

25.000 

.206803E+04 

26 

.206849E+01 

.694E-03 

24.000 

.186636E+04 

27 

.194847E+01 

.722E-03 

23.077 

. 151786E+04 

28 

.212907E+01 

.750E-03 

22.222 

. 193032E+0^ 

29 

.192044E+01 

.778E-03 

21.429 

.237375E+04 

30 

.181704E+01 

.806E-03 

20.690 

.157916E-^04 

31 

•  181628E-^01 

•833E-03 

20.000 

.89l445E-i-03 

TABLE  3 

THE  PGHHi  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIEU)  AT  ASHKHABAD 
20  SEPTEMBSt  1937  (LOW-PASS  FUTESED). 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

1 

.179439E+05 

•OOOE+OO 

.000 

2 

.292615E+04 

.lllE-03 

150.000 

3 

-.949157E+04 

.222E-03 

75.000 

4 

-.533074F+04 

.333E-03 

50.000 

5 

.101460E+04 

,444E-03 

37.500 

6 

.317897E+04 

.556E-03 

30.000 

7 

.171715E+04 

.667E-03 

25.000 

8 

-.960888E+02 

.778E-03 

21.429 

9 

-•144527E+04 

.889E-03 

18.750 

10 

-.150309E+04 

.lOOE-02 

16.667 

11 

.717671E+03 

.lllE-02 

15.000 

12 

.109009E+04 

.122E-02 

13.636 

13 

-.522785E+03 

.133E-02 

12.500 

14 

-.226249E+03 

.144E-02 

11.538 

15 

.133347E+04 

.156E-02 

10.714 

16 

.570191E+03 

.167E-02 

10.000 

17 

-.628134E+03 

.178E-02 

9.375 

18 

-.139616E+04 

.189E-02 

8.824 

19 

-.106033E+04 

.200E-02 

8.333 

20 

.194110E+03 

.211E-02 

7.895 

21 

.457800E+03 

.222E-02 

7.500 

22 

.113025E+04 

.233E-02 

7.143 

23 

.826872E+03 

.244E-02 

6.818 

24 

-.176093E+04 

.256E-02 

6.522 

25 

-.216693E+04 

•267E-02 

6.250 

26 

.537400E+03 

.278E-02 

6.000 

27 

.245529E+04 

.289E-02 

5.769 

28 

.676671P+03 

.300E-n2 

5.556 

29 

-.183471E+04 

.311E-02 

5.357 

30 

-.513085E+03 

.322E-02 

5.172 

31 

•377067E+02 

.333E-02 

5.000 

TABLE  4 

OSE  POWER  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIELD 
20  SEPTaiBER  19^7  (HIGH -PASS  FILTERED). 


POWER  DENSITY 


•OOOOOOE+OO 

.897712E+04 

•295033E+04 

•155729E+04 

.196698E+04 

.220064E+04 

.154380E+ft4 

.100167E+n4 

.100586E+04 

.742043E+03 

.426244E+03 

.319017E+03 

.254073E+03 

.179531E+03 

.131611E+03 

.106342E+03 

.112996E+03 

.106757E+03 

.788724E+02 

.574017E+02 

•417640E+02 

.442818E+02 

•515664E+02 

.438601E+02 

.363167E+02 

.376139E+n2 

.328361E+02 

.259808E+02 

.257291E+02 

.3n6650E+02 

.157647E+02 


AT  ASHKABAD 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

.249597E+02 

,000E+00 

.000 

.OOOOOOE+00 

2 

.200757E+02 

•132E-03 

126.400 

.223627E+05 

3 

•143485E+02 

•264E-03 

63.200 

.217677E+05 

4 

.117722E+02 

.396E-03 

42.133 

.104099E+05 

5 

.103865E+02 

.527E-03 

31.600 

.428982E+04 

6 

.810345E+01 

•659E-03 

25.280 

.327285E+04 

7 

.656234E+01 

.791E-03 

21.067 

.339714E+04 

8 

.550026E+01 

.923E-03 

18.057 

.400704E+04 

9 

.421846E+01 

.105E-02 

15.800 

.356261E+04 

10 

.204951E+01 

.119E-02 

14.044 

.211246E+04 

11 

.224725E+00 

.132E-02 

12.640 

.145312E+04 

12 

-.155644E+01 

.145E-02 

11.491 

.214154E+04 

13 

-•359934E+01 

.158E-02 

L0.533 

.243596E+04 

14 

-.543769E+01 

.171E-02 

9.723 

.149564E+04 

15 

-.541180E+01 

.185E-02 

9.029 

.120960E+04 

16 

-•403015E+01 

.198E-02 

8.427 

.218324E+04 

1  7 

-•462940E+01 

.211E-02 

7.900 

.233770t+04 

18 

-•632479E+01 

.224E-02 

7.435 

.135818E+04 

19 

-•672586E+01 

.237E-02 

7.022 

.825243E+03 

20 

-•599494E+01 

.251E-02 

6.653 

.941175E+03 

21 

-•565225E+01 

.264E-02 

6.320 

.876061E+03 

22 

-.549113E+01 

.277E-02 

6.019 

.528636E+03 

23 

-.513436E+01 

.290E-02 

5.745 

.429583E+03 

24 

-.435678E+01 

•303E-02 

5.496 

.543620E+03 

25 

-•362216E+01 

.316E-02 

5.267 

.548557E+03 

26 

-•315495E+01 

.330E-02 

5.056 

.405881E+03 

27 

-•277284E+01 

.343E-02 

4.862 

.307463E+03 

28 

-.218406E+01 

.356E-02 

4.681 

.472505E+03 

29 

-.201542E+01 

.369E-02 

4.514 

.878581E+03 

30 

-.224663E+01 

.382E-02 

4.359 

.116149E+04 

31 

-.174977E+01 

.396E-02 

4.213 

.655667E+03 

TABLE  5 

THE  PGHH)  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIELD  AT  LVOV 
20  SEFTESfBER  1997  (LOW-PASS  FILTESED). 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

.542315E+05 

•OOOE+OO 

.000 

.OOOOOOE+00 

2 

.222632E+04 

.527E-03 

31.600 

.529469E+04 

3 

-.302254E+05 

.105E-02 

15.800 

.252012E+04 

4 

-.831475E+04 

.158E-02 

10.533 

.178999E+04 

5 

.968791E+04 

•211E-02 

7.900 

.113665E-f04 

6 

.124626E+05 

.264E-02 

6.320 

.599837E+03 

7 

-.184197E+04 

.316E-02 

5.267 

.405788E+03 

8 

-.768729E+04 

.369E-02 

4.514 

.399254E+03 

9 

-.480041E+03 

.422E-02 

3.950 

.316888E+03 

10 

.186293E+04 

.475E-02 

3.511 

.207710E+03 

11 

-.202412E+04 

.527E-02 

3.160 

.184046E+03 

12 

,456927E+03 

.580E-02 

2.873 

.175130E+03 

13 

.286874E+03 

.633E-02 

2.633 

.137327E+03 

14 

-.123826E+04 

.686E-02 

2.431 

.128924E+03 

15 

-.162086E+03 

.738E-02 

2.257 

.135249E+03 

16 

-.199699E+04 

.■^91E-02 

2.107 

.102847E+03 

17 

.936901E+03 

•844E-02 

1.975 

.750166E+02 

18 

.531504E+04 

.P97E-02 

1.859 

.627144E+C2 

19 

.310242E+03 

.949E-02 

1.756 

.4522fl6E+02 

20 

-•638295E+04 

.lOOE-01 

1.663 

.370239E+02 

21 

-.276170E+04 

.105E-01 

1.580 

.458348E+C2 

22 

.527089E+04 

.1 llE-01 

1.505 

.561107E+02 

23 

.393865E+04 

.116E-01 

1.436 

.501630E+02 

24 

-.122666E+04 

.121E-01 

1.374 

.409907E+02 

25 

-.650960E+04 

.127E-01 

1.317 

.349222E+02 

26 

.123383E+03 

.132E-01 

1.264 

.346569E+02 

27 

.588585E+04 

.137E-01 

1.215 

.371423E+C2 

28 

-.887667E+03 

.142E-01 

1.170 

.294143E+02 

29 

-.212820E+04 

.148E-01 

1.129 

.2357e6E+02 

30 

-.221547E+04 

.153E-01 

1.090 

.188827E+02 

31 

•129949E+04 

•158E-01 

1.053 

.970753E+01 

TABLE  6 

THE  POWiH  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIELD  AT  LVOV 
20  SEPTPIBm  19^7  (HIGH-PASS  FILTraSD). 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

,527668E+00 

•OOOE+OO 

.000 

.OOOOOOE+00 

2 

•452818E400 

.881E-04 

189.200 

.976567E+03 

3 

•360401E+00 

.176E-03 

94.600 

.298734E+03 

4 

.308487E+00 

•264E-03 

63.067 

.119504E+03 

3 

•279587E+00 

.352E-03 

47.300 

.152760E+03 

6 

•249109E+00 

.440E-03 

37.840 

.117310E+03 

7 

•249223E+00 

.529E-03 

31.533 

.101320E+03 

8 

.266438E+00 

.617E-03 

27.029 

.137529E+03 

9 

.266475E+00 

.705E-03 

23,650 

.111521E+03 

10 

.246519E+00 

.793E-03 

21.022 

.545452E+02 

11 

,231296E+00 

•881E-03 

18.920 

.351943E+02 

12 

•  218261E-f00 

.969E-03 

17.200 

.415971E+02 

13 

.199516E+00 

•106E-02 

15.767 

.375442E+C2 

14 

•177892E+00 

•115E-02 

14.554 

.221671E+C2 

15 

.187680E+00 

.123E-02 

13.514 

.227775E+C2 

16 

•207658E+00 

.132E-02 

12.613 

.460398E+02 

17 

.197242E+00 

.141E-02 

11.825 

.575541E+02 

18 

.158906E+00 

•150E-02 

11.129 

.347276E+02 

19 

.127072E+00 

.159E-02 

10.511 

.186129E+02 

20 

•100968E+00 

.167E-02 

9.958 

.204121E+C2 

21 

.632944E-01 

.176E-02 

9.460 

.161251E+02 

22 

•423435E-01 

.185E-02 

9.010 

.909478E+01 

23 

.394748E-01 

.194E-02 

8.600 

.312655E+01 

24 

.386641E-01 

.203E-02 

8.226 

.779160E+01 

25 

.384960E-01 

.211E-02 

7.883 

.951384E+01 

26 

•481945E-01 

.220E-02 

7.568 

.141426E+02 

27 

.614213E-01 

•229E-02 

7.277 

.139312E+02 

28 

.626074E-01 

.238E-02 

7.007 

.114823E+02 

29 

.436504E-01 

•247E-02 

6.757 

.172622E+C2 

30 

.146695E-01 

•255E-02 

6.524 

.271249E+02 

31 

-.190410E-01 

.264E-02 

6.307 

.153122E+02 

TABLE  7 

THE  POWER  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  PIEID  AT  TBILESI 
20  SEPTEMBER  19^7  (LOW-PASS  FILTERED). 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

.521954E+03 

.OOOE+00 

•  000 

.OOOOOOE+00 

2 

.521579E+02 

.352E-03 

47.300 

.126653E+03 

3 

-,289376E+03 

.705E-03 

23.650 

.578031E+02 

4 

-.958594E+02 

.106E-02 

15.767 

.359559E+02 

*) 

.127578E+03 

.141E-02 

11.825 

.259795E+02 

6 

."998306E+02 

.176E-02 

9.460 

.140052E+02 

7 

-.195070E+02 

.211E-02 

7.883 

.831244E+al 

8 

-.597786E+b2 

.247E'-02 

6.757 

.775461E+01 

9 

-,494442E+02 

.282E-02 

5.913 

.456011E+01 

10 

-;5'8'80b6E+01 

.317E-02 

5.256 

.234601E+0i 

11 

.564351E+02 

.352E-02 

4.730 

.222932E+01 

12 

.211993E+02 

.388E-02 

4.300 

.277142E+01 

13 

-.492906E+02 

.423E-02 

3.942 

•260917E+01 

14 

-•496180E+b2 

.458E-02 

3.638 

.188579E+01 

15 

-.968609E+01 

.493E-02 

3.379 

•153600E+01 

16 

.475798E+02 

.529E-02 

3.153 

.140298E+01 

17 

.364598E+02 

.564E-02 

2.956 

.122181E+01 

18 

-.263554E+01 

.599E-02 

2.782 

.109477E+01 

19 

-•428643E+02 

.634E-02 

2.628 

.876305E+00 

20 

-.36bi34E+02 

.669E-02 

2.489 

.606516E+bb 

21 

.183575E+02 

.705E-02 

2.365 

.443000E+00 

22 

.281596E+02 

•740E-02 

2.252 

.366511E+00 

23 

-.302384E+02 

.775E-02 

2.150 

.334333E+00 

24 

-.320544E+02 

.810E-02 

2.057 

.424907E+00 

25 

-.110150E+02 

.846E-02 

1.971 

.437602E+00 

26 

'.330273E+02 

.381E-02 

1.892 

.347043E+bb 

27 

.206792E-^02 

.916E-02 

1.819 

.315703E+00 

28 

-.466457E+02 

.951E-n2 

1.752 

.348084E+00 

29 

-.353344E-^02 

.987E-02 

1.689 

.450148E+O0 

30 

.371368E+02 

.102E-01 

1.631 

.437743E+00 

31 

•563057E+02 

.106E-01 

1.577 

.225042E+00 

TABLE  8 

THE  pares  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIELD  AT  TBILISI 
20  SEFTQffiES  19^7  (HIGH-PASS  FILTERED). 
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densities  are  -6.5  dlD  at  Ashkhabad  and  -10,8  db  at  Tbilisi.  These  relative 
power  levels  reflect  the  differences  in  local  topography  at  each  position. 

Magnetotelluric  Field  at  Tbilisi 

The  records  made  at  Tbilisi  on  1  and  2  September  1957  were  subjected 
to  detailed  analysis.  The  results  are  presented  in  Tables  9  through  I8  and 
Figures  21  through  29*  In  particular,  the  power  spectra  of  the  North-South 
and  East -West  components  of  the  magnetic  and  telluric  fields  at  Tbilisi  are 
computed  for  the  two  ten-hour  periods;  0600  hours  to  160O  hours  (GMT)  on 
1  September  1957  and  O5OO  hours  to  15OO  hours  (GMT)  on  2  September  1957*  An 
additional  computation  was  made  for  the  East -West  telluric  component  for  the 
twenty-hour  period  0600  hours  (GMT)  on  1  September  1957  to  0200  hours  (GMT) 
on  2  September  1957*  A  comparison  between  the  power  spectra  for  the  twenty- 
hour  sample  (shown  in  Figure  29)  and  the  ten -hour  saunples  (shown  in  Figures  25 
and  2U)  does  not  reveal  significant  difference  between  the  ten-  and  twenty- 
ho\ir  samples. 

For  all  of  the  1  and  2  September  data,  power  spectra  were  computed  using 
the  original  sampling  interval  (72  seconds)  and  then  power  spectra  were  com¬ 
puted  after  low -pass  filtering  and  decimation  to  a  new  sampling  interval  of 
288  seconds.  For  some  reason  that  cannot  be  explained,  when  the  East -West 
component  of  r,he  magnetic  field  and  the  North-South  component  of  the  telluric 
field  were  analyzed  with  the  aid  of  a  low -pass  filter,  the  power  spectra  were 
erratic  and  ^anreiiabJ e.  Bus  showed  up  botti  in  numerous  negative  values  for 
the  power  spectral  densities  and  large  values  of  coherency  (see  Tables  15 
and  17).  Tliese  data  were  rejected. 

The  spectra  ol  the  North-South  magnetic  and  East-West  telluric  components 
do  not  show  any  striking  depart’ures  in  behavior  from  the  curves  for  20  Septem¬ 
ber.  The  empirical  formula  in  equation  (58)  can  again  be  made  to  fit  the 
magnetic  spectra  with  n  =  2  and  the  teilurlc  spectra  with  a  slope  of  n  between 
1.1  and  1.2. 

A  comparison  of  the  power  levels  at  J  cpks  shows  that  the  variation  in 
any  component  from  1  to  2  September  is  less  than  1  db. 
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FIGURE  21 

THE  POWER  SPECTRUM  OF  NORTH-SOUTH  MAGNETIC  FIELD 
AT  TBILISI  I  SEPTEMBER  1957 


K 


U)  m  ^  ro 

SdO/gX  Nl  A1ISN30  d3/IA0d  901 


-  6s  - 


-  64  - 


-  65  - 


in 


^  to  <vj 

SdO/^X  Nl  AilSNBQ  d3M0d  901 


-  66  - 


FIGURE  27 

THE  POWER  SPECTRUM  OF  NORTH-SOUTH  TELLURIC  FIELD 
AT  TBILISI  I  SEPTEMBER  1957 
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FIGURE  29 

THE  POWER  SPECTRUM  OF  EAST-WEST  TELLURIC  FIELD 
AT  TBILISI  1-2  SEPTEMBER  1957 


AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

•111927E+01 

•  OOOE-fOO 

.000 

.OOOOOOE+00 

2 

.964989E+00 

.579E-04 

286.000 

.306533E+04 

3 

.813778E+00 

.116E-03 

144.000 

.173019E+04 

4 

.706746E+00 

.174E-03 

96.000 

.739613E+03 

5 

•636353E+00 

•231E-03 

72.000 

.414087E+03 

6 

•557655E+00 

.289E-03 

57.600 

.353632E+03 

7 

•502514E+00 

.347E-03 

48.000 

.257352E+03 

8 

•460640E+00 

.405E-03 

41.143 

.195588E+03 

9 

•386045E+00 

.463E-03 

36.000 

.230166E+03 

10 

.347015E+00 

.521E-03 

32.000 

.171586E+03 

11 

.293465E+00 

.579E-03 

26.800 

.115192E+03 

12 

•235451E+00 

.637E-03 

26.182 

.121073E+03 

13 

.197398E+00 

.694E-03 

24.000 

.114374E+03 

14 

.167107E+00 

.752E-03 

22.154 

.101032E+03 

15 

•139452E+00 

.810E-03 

20.571 

.887191E+02 

16 

.966622E-01 

.868E-03 

19.200 

.815121E+02 

17 

•631499E-01 

.926E-03 

18.000 

.767610E+02 

16 

.917138E-02 

.984E-03 

16.941 

.869785E+02 

19 

-.221510E-01 

.104E-02 

16.000 

.976911E+02 

20 

-.321277E-01 

.llOE-02 

15.158 

.874828E+02 

21 

-•184531E-01 

.116E-02 

14.400 

.575089E+02 

22 

.648993E-02 

.122E-02 

13.714 

.303470E+02 

23 

•158395E-01 

•127E-02 

13.091 

.234011E+02 

24 

.123592E-01 

.133E-02 

12.522 

.302503E+02 

25 

-•179451E-01 

.139E-02 

12.000 

.435613E+02 

26 

-.608128E-01 

.145E-02 

11.520 

.825516E+02 

27 

-.966493E-01 

.150E-02 

11.077 

.122263E+03 

28 

-.110420E+00 

.156E-02 

10.667 

.109925E+03 

29 

-.784250E-01 

.162E-02 

10.286 

.328892E+02 

30 

-.721355E-01 

.168E-02 

9.931 

.887128E-t-02 

31 

-•612256E-01 

.174E-02 

9.600 

.500789E+02 

TABLE  9 

OSE  TOUm  SPECroUM  FOR  THE  EAST-WEST  TELLURIC  FIELD  AT  TBILISI 
1-2  SEPTQffiBt  1997  (LCW-PASS  FII/MtED). 
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AUTOCORRELATION 

FREQUENCY 

PERIOD 

POWER  DENSITY 

1 

.498126E+04 

.000E-»-00 

•  000 

.OOOOOOE+00 

2 

-.372095E+03 

.231E-03 

72.000 

.414487E+03 

3 

-.300428E+04 

.463E-03 

36.000 

.171882E+03 

4 

-.145405E+03 

.694E-03 

24.000 

.103381E+03 

5 

.165148E+04 

.926E-03 

18.000 

•  708166E-t-02 

6 

-.758578E+02 

•116E-02 

14.400 

.565561E+02 

7 

-.100168E+04 

.139E-02 

12.000 

.570442E+02 

8 

.938160E+03 

.162E-02 

10.286 

.647195E+02 

9 

-.980898E+02 

.185E-02 

9.000 

.574098E+02 

10 

-.297661E+03 

•208E-02 

8.000 

.376909E+02 

11 

-.215272E+02 

•231E-02 

7.200 

.270151E+02 

12 

.102245E-«-01 

•255E-02 

6.545 

.219163E+02 

13 

.189779E+02 

.278E-02 

6.000 

.165020E+02 

14 

-•.853661E+02 

.301E-02 

5.538 

.186011E+02 

15 

.609078E+03 

.324E-02 

5.143 

•  281894E'^02 

16 

-.882691E+03 

.347E-02 

4.800 

.266271E+02 

17 

-.443019E+02 

.370E-02 

4.500 

.236774E+02 

18 

•930281E+03 

.394E-02 

4.235 

.310236E-i'02 

19 

-•286435E+02 

.417E-02 

4.000 

•242912E402 

20 

-.104106E+04 

.440E-02 

3.789 

.988067E+01 

21 

-.129441E+03 

.463E-02 

3.600 

.627292E+01 

22 

.141566E+04 

•486E-02 

3.429 

.488245E-I-01 

23 

-.510356E+03 

,509E-02 

3.273 

.319039E+01 

24 

-.706014E+03 

.532E-02 

3.130 

.287962E+01 

25 

.363506E+03 

.556E-02 

3.000 

.562041E+01 

26 

.420156E+03 

.579E-02 

2.880 

.166903E+02 

27 

-.379849E+03 

.602E-02 

2.769 

.222770E+02 

28 

-.404489E+03 

.625E-02 

2.667 

.110764E+02 

29 

.740547E+03 

•648E-02 

2.571 

.389341E+01 

30 

-•229801E+03 

•671E-02 

2.483 

.558628E-t>0l 

31 

-.592975E+02 

•694E-02 

2.400 

•  287188E<«>01 

TABLE  10 

THE  POMIR  SPECTRUM  FOR  THE  EAST-WEST  TELLURIC  FIELD  AT  TBILISI 
1-2  SEPTEMBER  1957  (HIGH-PASS  FILTERED). 
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EAST  NORTH 

FREQUENCY  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 
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;ghs  for  the  east-west  mai 
ramat  19^7  (hicb-pass  fe 


NORTH  EAST 

FREOUENCr  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 
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,  AND  mSK  REIATIOHS  FOR  1SE  BCRIB-SOU! 
AT  TBUJSI  1  SSraffiBI  19^7  (HKa-PASS 
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EREOOENCY  MAGNETIC  TELLURIC  COHtRfcNCr  PHASE  ANGLE 

CPS  PO<(ER  DENSITY  POWER  DENSITY  RADIANS 
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NOMIti  CAST 

FREQUENCY  MAbNEIiC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  poncR  density  power  density  radians 


(0’<0 

eo»oot^mtfMn  tf>  ^  oe^F4m<^^ir>o<v^  io.«o 


^^ll^u^rsl^r^J(nr^ 


I 


I  I 


^w4t^ 

I  I 


•  I  I 


I  I  ^  #H  pH  fH 


s 


w 

'T  *t  1^1 

lA 

m  m 

m  m 

OT 

m 

m 

CA 

P'1 

m 

<*1 

M 

PI 

CM 

f>J 

nj 

w 

nj 

p> 

P> 

p> 

PI 

m 

P3 

O 

o  o  o 

o 

C5 

o  o 

o  o 

CJ 

f  > 

C'> 

'  > 

r  ) 

'  ) 

C) 

'  ) 

f  ■> 

<  > 

'  > 

L  » 

'  3 

' ) 

■  ) 

■  "I 

'  J 

' ) 

*  T 

' ) 

3 

♦ 

♦  ♦  + 

4- 

4 

4  4 

4  4 

N* 

4 

t- 

1" 

H 

»• 

I- 

!• 

»- 

1- 

»• 

!• 

•' 

(• 

h 

Ui 

UJ  UJ  UJ 

lU 

U 

U  U 

U  UJ 

IJ 

•  J 

U 

U 

'  J 

U 

;  J 

J 

'  J 

J 

J 

J 

J 

■  i 

■  J 

•  1 

•  J 

J 

J 

J 

J 

o 

h-  lA  «0 

m 

O' 

o  ^ 

00 

O' 

o 

r- 

rA 

r- 

o 

in 

CO 

CO 

.n 

rg 

•H 

■'^ 

O 

03 

o 

rW 

r  • 

( 

J' 

o 

O  a>  tA 

^  m 

m  *r 

<M 

rg 

r  J 

•O' 

eo 

r- 

O' 

A 

O' 

13 

CD 

a' 

g3 

03 

o 

o 

«0  O  h- 

»r 

«  pH 

O 

> 

•T 

M 

■t 

o 

0) 

pH 

o 

pH 

-H 

33 

IJ 

r*', 

O 

r* 

o 

«A  <0  O) 

m 

rj  (A 

00 

o 

r-- 

tff 

o 

in 

g) 

00 

O' 

0- 

x> 

-H 

•  ‘ 

;a 

X. 

P'V 

'■g 

■'J 

o 

m  ^  « 

'C 

» «e 

O' 

O' 

-H 

pN 

iH 

♦ 

lA 

M 

CD 

ni 

p» 

'X 

m 

m 

« 

c- 

<A  IM  « 

'T 

Ai 

-H 

rH 

-H 

pH 

pH 

-H 

■  H 

« 

XI 

•r 

A 

fH 

pg 

rg 

-H 

•H 

'V 

§ 


o 

^  ^  o 

M'l 

‘t 

n 

n 

n 

t 

t 

•A 

A 

A 

A 

t 

4 

t 

A 

4  4 

4 

4 

4 

'A 

4  4 

lA 

A 

A 

o 

o  o  o 

o 

o 

o 

o 

3 

o 

D 

o 

O 

O 

O 

“3 

■3 

O 

o 

3 

O 

J 

Y 

; 

3 

3  1; 

o 

3 

4 

4  4  f 

f 

4 

p 

P 

f 

t- 

♦■ 

t 

»• 

r 

t* 

f- 

r 

•  r 

t* 

h 

I* 

1 

h  1- 

• 

► 

UJ 

Ui  UJ  tu 

‘U 

J 

'  J 

a 

J 

J 

J 

J 

J 

J 

U 

■J 

J 

u 

j 

J 

J  J 

J 

1 

ij  -J 

J 

J 

J 

o 

pH  go  ^ 

O' 

00 

A 

fM 

A 

■X 

pH 

-H 

T 

■  H 

vO 

"H 

O 

4- 

A 

4 

.r  x 

O 

70 

O 

4 

A 

o 

X 

0 

o 

pH  O  O 

4* 

p> 

pH 

ng 

lA 

(A 

eg 

rg 

!A 

fA 

lA 

4- 

O' 

A  rg 

rg 

A 

*  H 

pH  4*  O 

4 

X 

lA 

o 

pg  O  CO 

•X 

O' 

P> 

pg 

rg 

o 

0 

'A 

g 

«H 

»A 

■a 

O 

A 

O  pH 

> 

A 

4 

*■* 

'4  'H 

X 

r- 

pH 

o 

r*  ^  « 

lA 

fh 

pH 

fA 

« 

■p 

M 

O 

A 

A 

JD 

4- 

:A 

g 

-H 

>  -H 

N 

o 

D 

O 

>  -H 

o 

4 

■“p 

o 

O'  O  M 

♦V 

lA 

lA 

O' 

<0 

«A 

o 

o 

Al 

O 

(A 

4- 

O  CD 

«H 

Ai 

pH  O  O' 

o 

A 

o 

o 

ng  lA  pH 

CA 

O'  -H 

•H 

•H 

O' 

pH 

pH 

*  « 

pH 

ng 

pg 

pH 

'  H 

^4  .  H 

•  4 

4 

X 

rg  4 

pH 

^g 

4 

I  I 


O 

4 

PA 

CA 

fA 

<A 

ft 

<A 

3 

3 

tA 

lA 

o 

PA 

4A 

lA 

(A 

tA 

rg 

eg  <g 

<  g 

tg 

<g 

f  g 

rg 

(g 

rg 

rg 

rg 

(g 

o 

o 

o 

O 

C> 

O 

o 

O 

*3  *3 

O 

■> 

*  ) 

^  ) 

3 

'  3 

.  3 

3 

' 

5 

M 

UJ 

UJ 

UJ 

UJ 

UJ 

LJ 

J 

•J  J 

J 

-1 

J 

t  J 

?  J 

J 

J 

J 

J  J 

J 

•  J 

_ 

J 

o 

p«. 

Aw 

>0 

A 

4 

rg 

pH 

O  rO 

r*- 

-o 

4 

rO 

iM 

•H 

O' 

X 

'g 

O  r- 

.A 

A 

pH 

,T 

r- 

wA 

rg 

■D 

rO 

»0 

o 

tt 

PA 

pH 

O' 

A- 

A 

lA  O 

a 

o 

4 

rg 

O 

A 

fA 

-t 

O  -A 

pH 

'p 

'A 

•V 

4 

O 

O 

rg 

^A 

o 

*A 

X 

4 

O 

'0  rg 

*■» 

A 

O' 

.n 

pH 

o 

rg 

'jO 

o 

pH  pH 

rg 

‘  g 

.3 

'A 

4 

A 

♦  3 

O 

o 

lA 

pH 

rg 

rg 

A 

4 

4  A 

A 

sO 

■o 

■* 

a> 

0 

o 

pH 

•  H  -H 

p  -4 

■  4 

.  1 

.  1 

.  1 

.  -4 

f  1 

1  H 

pH 

•  4 

o^cqo^Oph 


'n  ^  ^  4)  flO  <>  a3  r^  •<^  •♦  A 


eg  ^  *g 


O'  O  pH 
CM  A 


75 


AT  TBIIJSI  1  SSTBIBHt  1957  ( 


EAST  NORTH 

FREQUENCY  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 


Owoo^p-<ojh>«rM0^4’^h-^ao<Hr-%or^^«-4aor^r^aoo€0«o(Mr^(^ 

Of«>aOf^^tf>0^tf\Or^«OCO^^OOr^ 


I  I  ^  ^  ^  I 

I  I  I 


I  1^ 


I 


Or4flDmmh*>^mf>j^fnfl0^wr^TOO^r^tf>4’«or^«D<ncoo^Ov  ^  otto 


o 

<M 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

0-4 

0-4 

0-4 

0-4 

1-4 

o 

o 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

o 

♦ 

+ 

+ 

♦ 

+ 

+ 

+ 

+ 

♦ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

4 

•f 

+ 

+ 

♦ 

•f 

♦ 

+ 

+ 

+ 

+ 

IIJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

Ul 

Ul 

UJ 

UJ 

UJ 

Ul 

c 

a> 

O' 

f*> 

yO 

O' 

iTk 

«l> 

'O 

O' 

CM 

<^ 

CM 

tn 

00 

'O 

O' 

CM 

O' 

yC 

CO 

CM 

o 

rH 

CO 

'0 

f- 

o 

CM 

o 

o 

CM 

O' 

r- 

•-< 

O' 

O' 

lO 

CO 

'O 

CM 

CO 

'O 

CO 

'O 

r'- 

CO 

in 

O 

00 

CO 

CO 

CO 

yO 

o 

CM 

00 

K 

CM 

4- 

'O 

lA 

o 

CM 

CM 

CM 

O' 

O' 

e'¬ 

CM 

<0 

h- 

'O 

CM 

CM 

CM 

o 

CM 

CO 

r- 

O' 

O' 

m 

yO- 

o 

f^ 

m 

CM 

lO 

ir\ 

CM 

<^ 

en 

CM 

CO 

0-4 

o 

in 

r-4 

'O 

o 

oo 

0-4 

CO 

Oy 

o 

to 

o 

O' 

ii> 

CO 

OD 

'O 

CM 

CM 

m 

CM 

O' 

'0 

r- 

CO 

0-4 

o 

yO 

o 

r- 

o 

• 

00 

• 

• 

CO 

• 

CM 

• 

CM 

• 

CM 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CM 

• 

• 

f- 

• 

00 

• 

m 

• 

00 

• 

h- 

• 

r- 

• 

m 

• 

• 

• 

• 

yO 

• 

CO 

• 

co 

• 

» 


o 

m 

in 

4* 

CO 

<o 

<o 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

O 

O 

o 

O 

O 

O 

O 

O 

O 

o 

O 

o 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

UJ 

Ul 

Ul 

Ul 

Ul 

UJ 

Ul 

Ul 

UJ 

UJ 

Ul 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

UJ 

Ul 

UJ 

UJ 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

UJ 

o 

o 

00 

00 

CM 

in 

c^ 

m 

in 

CM 

O' 

CM 

00 

rH 

«-4 

4* 

o 

O' 

'O 

yC 

O' 

in 

CM 

O' 

m 

o 

o 

'O 

00 

o 

>0 

CM 

O 

CO 

CM 

CO 

CM 

CO 

CO 

O' 

CO 

O 

0-4 

in 

O' 

0-4 

m 

f^ 

CO 

CM 

o 

0-4 

o 

CM 

CO 

M) 

00 

o 

CM 

r- 

0-4 

O' 

m 

CM 

h* 

'O 

m 

00 

00 

CM 

w4 

CO 

m 

<• 

O' 

n- 

n- 

CO 

CO 

m 

O' 

yO 

o 

CM 

yO 

O' 

0-4 

yO 

CM 

O' 

oo 

yC 

in 

<o 

'0 

O' 

yt 

CO 

4- 

CO 

CM 

^4 

OD 

4- 

yO 

00 

0-4 

Oy 

o 

in 

m 

y* 

CM 

'O 

O' 

CM 

CO 

O 

in 

00 

m 

00 

00 

CM 

CM 

CM 

O 

fO 

O' 

Qy 

o 

O' 

yO 

CO 

o 

• 

m 

• 

• 

in 

• 

CM 

• 

0-4 

• 

0-4 

• 

r- 

• 

h- 

• 

O' 

• 

n- 

• 

-o 

• 

CO 

• 

CO 

• 

CM 

• 

• 

• 

• 

• 

O' 

• 

00 

• 

00 

• 

yC 

• 

• 

• 

CO 

• 

CO 

• 

m 

• 

in 

• 

in 

• 

m 

• 

O  CO 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

o  o 

4  1 

o 

1 

o 

1 

o 

1 

O 

O 

O 

1 

O 

O 

O 

I 

O 

O 

1 

O 

O 

1 

O 

O 

O 

1 

O 

O 

O 

1 

o 

O 

1 

O 

1 

O 

O 

1 

o 

O 

O 

1 

O 

O 

1 

Ul  Ul 

Ul 

Ul 

1 

Ul 

Ul 

UJ 

1 

UJ 

UJ 

UJ 

1 

Ul 

UJ 

1 

UJ 

Ul 

1 

Ul 

UJ 

Ul 

1 

Ul 

Ul 

UJ 

1 

Ul 

UJ 

1 

UJ 

1 

Ul 

Ul 

1 

Ul 

UJ 

1 

Ul 

1 

UJ 

UJ 

1 

Ul 

o  m 

O 

4“ 

O' 

O' 

o 

CM 

CO 

in 

'O 

oo 

O' 

CM 

4 

in 

r- 

00 

o 

CO 

4 

>0 

n* 

O' 

o 

CO 

4 

o  ^ 

CO 

4 

m 

in 

OO 

CM 

m 

00 

4 

O 

4 

O 

CO 

yO 

O' 

CO 

yC 

O' 

CM 

m 

00 

in 

00 

0-4 

4 

O  CO 

'0 

O' 

CM 

CO 

yC 

00 

o 

CO 

m 

o 

CM 

4 

O' 

0-4 

CO 

yO 

00 

o 

CO 

m 

o 

CM 

4 

0- 

O' 

O  CM 

4- 

'O 

O' 

0-4 

0-4 

0-4 

CM 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

CO 

4 

4 

4 

4 

in 

in 

m 

m 

yO 

yO 

'0 

'O 

yO 

<^rym^iA^^«00«0«-N<Mc^<^ii>«or^«>0'0(HfM(v^<4‘tf>^r>coo^O«-4 


76 


IE  FOHHI  SPECISUM,  COBHiENCY,  AHD  PHASE  REIATIQNS  FOR  OHE  EAST-WEST 
HGBTH-SOUIH  TEIIURIC  FIELDS  AT  TBIUSI  2  SEFTQffiSi  1997  (HKS-PASS 


NORTH  EAST 

FREQUENCY  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 


oaiOOO^*<o^o^r^eo«o^•^NON^^•f'>^i^F4^lr>o«^o^09« 

||p4pi4|||  || 

I  I  I  I  1  II 


OO^rtf^K^4’N(^«or^J^>^iC^intf^^•^aooc^l4‘«oaoa>(<^rg4'^ao^m 

ONmmtf>r^NO*^r>r*>o^^ao4'09^o<^ojo^NO^^o^p-4Nr^^ 


o 

4> 

ro 

ff\ 

ri 

m 

fA 

(A 

AJ 

AJ 

CA 

fA 

*A 

AJ 

AJ 

<A 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

Al 

AJ 

Al 

AJ 

o 

O 

O 

o 

o 

o 

o 

o 

O 

o 

O 

o 

O 

O 

O 

o 

o 

o 

O 

o 

o 

O 

O 

O 

o 

o 

o 

O 

O 

O 

O 

■f 

+ 

•♦> 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

•f 

+ 

*f 

+ 

+ 

♦ 

(AJ 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

ux 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

uu 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ui 

UJ 

o 

«M 

>o 

00 

pH 

0^ 

00 

CM 

r* 

IT 

oo 

<A 

o 

lA 

>0 

lA 

00 

o 

(A 

«A 

o 

pH 

O' 

o 

O' 

00 

fSI 

N 

00 

(A 

00 

AJ 

AJ 

0^ 

o 

o 

O' 

O' 

lA 

AJ 

O' 

PH 

00 

«0 

00 

O 

4- 

r* 

pH 

o 

>c 

-t 

00 

>0 

0^ 

fH 

00 

AJ 

A- 

lA 

m 

<0 

-o 

OO 

O' 

'0 

(A 

00 

fA 

AJ 

4 

O' 

«A 

o 

oo 

o 

^• 

lA 

flO 

o 

o 

AJ 

o 

AJ 

o 

AJ 

(A 

fA 

lA 

lA 

P4 

-0 

00 

CA 

Al 

00 

o 

ra 

h- 

o 

00 

O' 

o 

AJ 

lA 

AJ 

lA 

Al 

AJ 

o 

pH 

fA 

A- 

n 

o 

00 

AJ 

fA 

O' 

lA 

lA 

o 

• 

• 

• 

• 

• 

• 

• 

N 

• 

• 

• 

o« 

• 

<0 

• 

• 

• 

h4 

• 

00 

• 

O' 

• 

• 

OO 

• 

• 

CA 

• 

CA 

• 

• 

• 

fA 

• 

Al 

• 

CA 

• 

fA 

• 

Al 

• 

fA 

• 

4 

• 

O 

4 

lA 

(A 

lA 

lA 

lA 

4 

4 

4 

4 

4 

4 

4 

4 

{A 

fA 

4 

CA 

fA 

fA 

(A 

(A 

fA 

fA 

fA 

fA 

(A 

fA 

fA 

fA 

O 

O 

o 

O 

O 

O 

o 

O 

O 

O 

O 

O 

O 

O 

O 

o 

o 

O 

o 

o 

o 

O 

O 

o 

o 

o 

o 

o 

O 

O 

o 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

Ul 

Ul 

UJ 

UJ 

UJ 

Ul 

UJ 

UJ 

UJ 

Ul 

Ul 

Ul 

Ul 

Ul 

Ul 

UJ 

Ul 

UJ 

Ul 

Ul 

Ul 

Ul 

o 

(A 

K 

O' 

pH 

00 

'O 

AJ 

o 

AJ 

AJ 

o 

(A 

'O 

(A 

o 

'0 

o 

4 

>0 

O' 

4 

Al 

Al 

pH 

o 

fA 

00 

4 

o 

O' 

o 

lA 

m 

o 

o 

O' 

O' 

o 

'O 

4A 

pH 

AJ 

A- 

(A 

pH 

pH 

AJ 

pH 

AJ 

pH 

fA 

O 

O' 

o^ 

O' 

fA 

00 

1^ 

o 

pH 

(A 

pH 

AJ 

pH 

4 

4 

lA 

AJ 

o 

4 

o 

pH 

(A 

4 

00 

4 

fA 

fA 

AJ 

00 

tA 

'O 

CA 

AJ 

O' 

o 

«A 

'O 

fA 

lA 

lA 

A 

lA 

4 

AJ 

lA 

o 

o 

O 

A* 

Al 

O 

O' 

O 

4 

O 

00 

OO 

4 

fA 

pH 

00 

00 

f'“ 

o 

lA 

00 

4 

pH 

pH 

o 

O' 

O' 

Al 

tA 

lA 

pH 

lA 

»A 

o 

O 

4 

O 

o 

O' 

O' 

AJ 

O' 

00 

o 

00 

lA 

o 

• 

pH 

• 

AJ 

• 

pH 

• 

pH 

• 

pH 

• 

PH 

• 

4 

• 

(A 

• 

AJ 

• 

pH 

• 

pH 

• 

pH 

• 

pH 

• 

pH 

• 

• 

00 

• 

• 

4 

• 

Aj 

• 

AJ 

• 

Aj 

• 

pH 

• 

• 

pH 

• 

Al 

• 

pH 

• 

pH 

• 

AJ 

• 

• 

pH 

• 

o 

CA 

fA 

fA 

fA 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

AJ 

Al 

Al 

AJ 

AJ 

AJ 

Al 

A 

A 

A 

A 

AJ 

A 

A 

A 

A 

A 

A 

o 

+ 

UJ 

o 

1 

O 

1 

O 

1 

O 

1 

o 

1 

O 

1 

o 

1 

O 

1 

O 

1 

o 

1 

O 

1 

O 

O 

1 

O 

1 

O 

1 

O 

1 

O 

o 

O 

4 

O 

O 

O 

1 

O 

o 

1 

O 

O 

O 

1 

O 

O 

o 

1 

Ul 

1 

Ui 

Ul 

1 

Ul 

1 

Ul 

U) 

1 

UJ 

1 

UJ 

1 

Ul 

1 

Ul 

1 

UJ 

Ul 

1 

Ul 

1 

UJ 

1 

UJ 

1 

Ul 

Ul 

UJ 

I 

Ul 

UJ 

Ul 

1 

Ui 

UJ 

1 

Ul 

UJ 

1 

Ul 

1 

UJ 

1 

Ui 

1 

Ul 

1 

UJ 

o 

lA 

o 

4 

O' 

A- 

O' 

o 

AJ 

fA 

»A 

'0 

oo 

O' 

pH 

Al 

4 

tA 

r- 

00 

O 

pH 

CA 

4 

'O 

r- 

O' 

o 

pH 

(A 

4 

o 

pH 

fA 

4 

lA 

«A 

00 

AJ 

lA 

OO 

pH 

4 

o 

4 

1^ 

O 

fA 

'O 

O' 

fA 

'O 

O' 

AJ 

tA 

00 

pH 

lA 

00 

pH 

4 

o 

fA 

O' 

Al 

pH 

fA 

'O 

00 

c 

<A 

lA 

o 

AJ 

4 

r- 

O' 

pH 

CA 

'O 

00 

o 

fA 

lA 

h* 

o 

A 

4 

O' 

o 

AJ 

4 

O' 

pH 

pH 

pH 

AJ 

AJ 

AJ 

AJ 

fA 

fA 

fA 

CA 

fA 

4 

4 

4 

4 

tA 

tA 

tA 

lA 

4 

'O 

'O 

'O 

O  1-4 
«0  ft\ 


77 


EAST  NORTH 

FREQUENCY  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 


H 


O<-<4‘^4'rHrg^O^O4’«>0O«04’'Hr^a0O9'«04’fl0r^  Of^^O^OO' 

I  I 


o  m  4-  flo 
o  ^  o  n 

O  O  0  O 

•  •  •  •  • 

CM  m  0 


O  if>  4r\  in  tTT 

o  o  o  o  o 

+  ♦♦  +  + 

UJ  UJ  Ui  UJ  LU 

C  rvj  nj  4  cn 

O  4  cn  cn  fo 

o  m  o  4 

O  4  CM  r-  fM 

o  pH  CO  <n 

O  «H  4  »H  CO 


I  II  I  I  I  I  I  I  I  I  I  III 


o 

r- 

1^ 

>0 

0 

m 

m 

in 

in 

m 

m 

in 

in 

m 

in 

*n 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

4 

in 

0 

0 

c 

o 

O 

a 

o 

o 

o 

o 

o 

O 

o 

O 

o 

o 

o 

c 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

o 

o 

O 

+ 

+ 

■f 

+ 

+ 

■f 

+ 

+ 

+ 

+ 

+ 

•f 

+ 

+ 

+ 

+ 

+ 

•f 

+ 

-f 

+ 

+ 

+ 

+ 

+ 

+ 

•f 

+ 

-f 

+ 

+ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

o 

cn 

m 

o 

fH 

0 

4 

in 

m 

m 

'O 

m 

4 

•M 

o 

O' 

<n 

in 

O' 

in 

4 

o 

tn 

0 

0 

O' 

cn 

<n 

0 

o 

O' 

4 

r- 

O' 

n- 

0 

0 

cn 

o 

r- 

>o 

O 

sO 

nj 

n- 

>o 

0 

0 

0 

nj 

O 

4 

r- 

O' 

0 

O' 

4 

0 

O' 

r- 

o 

0 

4 

m 

o 

«M 

m 

r*» 

O' 

O' 

r- 

m 

<n 

4 

4 

4 

O 

C^ 

O' 

O' 

0 

0 

rg 

0 

CO 

o 

0 

m 

0 

o 

o 

r* 

0 

0 

o 

O 

fO 

fVI 

rsi 

4 

O' 

O' 

4 

in 

4 

m 

O 

4i 

0 

in 

^- 

0 

nj 

0 

0 

o 

r- 

in 

O' 

0 

O' 

o 

4 

0 

ni 

ni 

nj 

0 

m 

O' 

in 

0 

o 

4 

O' 

4 

O 

4 

0 

o 

in 

fO 

0 

O' 

o 

• 

0 

• 

*-H 

• 

# 

• 

• 

m 

• 

Ki 

• 

• 

• 

• 

• 

• 

• 

(n 

• 

rsj 

• 

>0 

• 

rg 

• 

• 

• 

• 

in 

• 

0 

• 

4 

• 

4 

• 

4 

• 

4 

• 

0 

• 

w 

• 

4 

• 

0 

• 

I  I  I 


o 

4 

0 

0 

0 

0 

0  0 

0 

0 

0 

0 

0 

0 

0  0 

0 

0 

eg 

eg 

rg 

eg 

eg 

CNJ 

eg 

eg 

eg 

eg 

eg 

fSJ 

eg 

o 

o 

1 

o 

O 

1 

O 

o 

1 

o  o 

1  1 

O 

1 

O 

o 

1 

o 

1 

O 

o 

1 

o  o 

O 

1 

o 

1 

o 

o 

O 

O 

O 

1 

o 

O 

1 

O 

O 

1 

O 

O 

o 

O 

UJ 

1 

UJ 

UJ 

1 

UJ 

UJ 

1 

UJ 

1  1 
UJ  UJ 

1 

UJ 

UJ 

1 

UJ 

1 

UJ 

UJ 

1 

UJ 

U(  UJ 

1 

UJ 

1 

UJ 

Ul 

UJ 

UJ 

UJ 

UJ 

UJ 

1 

UJ 

Ut 

UJ 

Ul 

UJ 

UJ 

UJ 

o 

r- 

1^ 

0 

0 

4 

nj  -H 

o 

0 

e^ 

0 

4 

0 

eg  pn 

O' 

0 

eg 

o 

r«* 

0 

0 

pH 

O' 

0 

eg 

o 

0 

0 

o 

0 

0 

0 

O' 

f'*  0 

0 

o 

0 

0 

4 

eg 

O  0 

0 

0 

4 

o 

0 

pH 

r* 

0 

0 

4 

o 

0 

eg 

e- 

0 

o 

pH 

^- 

0 

0 

4  O 

0 

eg 

r- 

0 

O' 

0 

^  0 

eg 

0 

O 

pH 

eg 

eg 

0 

0 

4 

0 

0 

0 

0 

n- 

c 

• 

0 

• 

pH 

• 

pH 

• 

eg 

• 

eg 

• 

0  4 

•  • 

4 

• 

0 

• 

0 

• 

0 

• 

0 

• 

e* 

• 

0  0 

•  • 

O' 

• 

O' 

• 

• 

pH 

• 

pH 

• 

• 

pH 

• 

pH 

• 

pH 

• 

• 

pH 

• 

pH 

• 

pH 

• 

pH 

• 

pH 

• 

O' 

0 

o 

0  0 

pH 

O  0 

0 

r*  h- 

o 

O' 

O' 

0 

o 

0 

e- 

0 

O 

4 

O' 

0 

C 

4 

0 

4 

eg  0 

O 

0  0 

0 

0  0 

0 

O' 

o 

eg 

e« 

e- 

0 

o 

0 

o 

C 

0 

O' 

0 

eg 

0 

O  O 

0 

O'  0 

0 

eg  o 

fM 

4 

0 

O' 

eg 

O' 

o 

0 

4 

0 

• 

• 

• 

•  • 

• 

•  • 

• 

•  • 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

f- 

eg 

0 

4  ^ 

0 

0  eg 

eg 

2 

3 

r- 

o 

pH 

O 

O' 

0 

0 

pH 

0 

0 

m 

0 

0 

ni 

pH 

0 

pH 

pH 

•H 

rg 


0 

0 

4 

4  eg 

0 

0  0 

0 

0  0 

4 

0 

0 

4 

eg 

4 

4 

0 

4 

0 

0 

4 

0 

4 

4 

o 

o 

O 

O  O 

o 

o  o 

o 

o  o 

O 

o 

o 

O 

O 

o 

O 

o 

O 

O 

o 

O 

o 

O 

O 

+ 

+ 

+ 

+  + 

+ 

+  + 

+ 

+  + 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

UJ 

UJ 

UJ 

UJ  UJ 

UJ 

UJ  UJ 

UJ 

UJ  UJ 

UJ 

UJ 

UJ 

ui 

LU 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

O' 

0 

4 

0  eg 

eg 

0  pH 

eg 

O'  f^ 

0 

O' 

0 

0 

>o 

0 

0 

eg 

0 

0 

0 

o 

eg 

o 

0 

eg 

0 

0 

0  e- 

O  0 

O' 

0  eg 

0 

0 

O' 

O' 

0 

o 

eg 

0 

eg 

0 

0 

c 

0 

0 

r- 

eg 

O 

0  O 

0 

0  pH 

4 

0  e- 

h* 

4 

0 

0 

4 

4 

rg 

0 

4 

o 

4 

eg 

eg 

0 

4 

4 

eg 

0 

4  O 

4 

O'  o 

0 

4  0 

0 

PH 

4 

o 

0 

0 

0 

0 

eg 

O' 

4 

pH 

o 

0 

0 

4 

0  4 

pH 

O'  O' 

eg 

pH  CO 

0 

0 

0 

0 

e* 

4 

O' 

0 

0 

0 

0 

O 

4 

0 

pH 

eg 

0  0 

eg 

PH  0 

-H 

eg  pH 

pH 

pH 

pH 

0 

nj 

0 

4 

pH 

0 

pH 

O' 

I 

§ 

CO 


CO 

a 


•Hrgcr>4in^r^eo<>O^rMm4tn0r^«oo*»OpHnj<n4Xt0r*»0O'O^ 


78 


NORTH  EAST 

FREQUENCY  MAGNETIC  TELLURIC  COHERENCY  PHASE  ANGLE 

CPS  POWER  DENSITY  POWER  DENSITY  RADIANS 


fn 

CO 

h- 

tr> 

ID 

if>  O' 

O 

o 

4 

fM 

4 

O' 

kD 

fO 

4 

<0 

flO 

ir\ 

m 

m 

Ifl  »H 

CO 

CO 

4 

CO 

4 

4 

o 

CO 

o 

•H 

4 

fM  m 

tn 

4 

4  rg 

m 

O' 

IM 

CO 

/> 

4 

rg 

O' 

4 

4 

I  ^  ^  ^  ^  ^ 

III  I 


I  I  I  l-i^ 

I  III 


00'9«o^ocM^aocamOtOrt^^^r^«rM«<«o^m^rnco9«0(n 
o«oco9k0^o«^^OirtN>«^ir\^  «tf\oo«o<^^Choo«ooo9J>^rMm 


O 

iD 

4 

4 

4 

4 

CD 

(D 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

CD 

o 

o 

o 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

O 

o 

o 

o 

O 

O 

O 

O 

o 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

-f 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

lU 

uu 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

UJ 

LU 

UJ 

UJ 

UJ 

UJ 

OJ 

Ui 

UJ 

UJ 

UI 

UI 

UI 

UJ 

ul 

UI 

UJ 

UI 

UI 

UJ 

UI 

UI 

UI 

o 

fSI 

O' 

ID 

4 

4 

oo 

O' 

4 

CD 

r- 

iD 

^4 

r- 

4 

r- 

00 

t-4 

•sO 

O' 

CM 

p«- 

O' 

CD 

4 

c 

O' 

rH 

o 

4 

O 

CM 

o 

4 

O' 

CM 

CM 

4 

00 

CM 

O' 

O' 

CO 

r^ 

4 

o 

4 

4 

O' 

kD 

CD 

O' 

4 

CM 

•D 

o 

m 

CO 

<D 

o 

t-i 

4D 

o 

if\ 

O' 

o 

CM 

^- 

r- 

CD 

O' 

QO 

r'- 

4 

•H 

CD 

CM 

O 

O' 

4 

CO 

4 

o 

CM 

CD 

CO 

CD 

O' 

O' 

4 

4 

CD 

CO 

O' 

O' 

4 

f^ 

4 

CO 

4 

•H 

CD 

CD 

ID 

4 

fD 

4 

4 

4 

4 

o 

iD 

o 

4 

o 

O 

CO 

«n 

» 

CM 

CO 

ID 

O' 

4 

O 

CD 

O' 

CD 

CO 

o 

CM 

CM 

CO 

o 

kD 

O' 

o 

•H 

4 

CM 

-H 

CO 

4 

4 

4 

4 

4 

4 

4 

ID 

CD 

m 

CM 

CM 

CM 

CM 

CD 

4 

4 

4 

4 

iD 

4 

4 

4 

9 


O 

r- 

4 

4 

O 

4 

4 

4 

14 

4 

tD 

4 

4 

4 

4 

4 

4 

4 

■D 

4 

4 

4 

4 

4 

4 

4 

4 

4 

n 

4 

o 

o 

O 

O 

O 

O 

o 

o 

O 

o 

O 

O 

O 

o 

C 

O 

C 

o 

O 

o 

O 

O 

O 

o 

C 

O 

C 

O 

O 

o 

o 

+ 

•f 

+ 

+ 

+ 

+ 

+ 

4- 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

UJ 

UJ 

UI 

UI 

ul 

UJ 

ul 

Ul 

xl 

Ul 

Ul 

Ul 

Ul 

UJ 

ul 

Ul 

Ul 

Ul 

Ul 

UJ 

ul 

Ul 

OJ 

<iJ 

ul 

Ul 

Ul 

Ul 

‘Jj 

Ul 

ul 

o 

r- 

CM 

O' 

r- 

4 

00 

4 

M 

4 

4 

CM 

4 

r>^ 

X 

4 

CM 

CM 

o 

-H 

CM 

4 

o 

X 

X 

CM 

CM 

O' 

4 

CD 

o 

O' 

I'- 

CD 

CM 

CO 

4- 

(D 

4 

4 

f^ 

r«- 

0- 

4 

X 

r- 

CM 

CD 

> 

4 

o 

X 

O 

4 

4 

kT) 

O' 

X 

-H 

o 

o 

4 

CM 

CM 

4 

O 

CM 

4 

CD 

4 

CD 

4 

h- 

4 

O' 

O' 

X 

»— 1 

CD 

CM 

CD 

4 

4 

4 

o 

P~ 

O' 

O 

4 

4 

O 

CM 

O' 

CM 

CO 

X 

X 

X 

CM 

o 

f'- 

1— I 

X 

CD 

o 

r- 

CM 

X 

4 

r- 

o 

ID 

4 

4 

■D 

O 

O' 

CM 

CM 

CO 

O' 

p^ 

X 

CM 

PE, 

o 

4 

—w 

O' 

CM 

p*- 

X 

CM 

o 

4 

•M 

CM 

r--4 

CM 

(r\ 

CM 

(M 

4 

4 

-H 

X 

4 

TD 

X 

'D 

X 

(D 

4 

O 

4 

CD 

CD 

CD 

CD 

m 

m 

'D 

(D 

D 

<D 

rr\ 

CD 

<D 

<D 

CD 

CD 

CM 

CM 

CM 

CM 

CM 

M 

CM 

CM 

CM 

CM 

CM 

CM 

M 

o 

o 

1 

o 

1 

O 

1 

O 

1 

O 

1 

o 

1 

o 

1 

O 

1 

O 

O 

1 

O 

O 

O 

O 

o 

O 

O 

O 

O 

o 

O 

O 

O 

O 

■Ti 

O 

O 

O 

O 

O 

UJ 

UJ 

Ul 

1 

Ul 

UJ 

1 

Ul 

1 

UJ 

1 

Ul 

uJ 

Ul 

1 

Ul 

UJ 

UJ 

u 

Ul 

Ul 

UJ 

aJ 

Ul 

UJ 

uJ 

UJ 

ul 

UJ 

ul 

Ui 

Ul 

Ul 

Ul 

Ul 

Ul 

o 

1^ 

r- 

4 

4 

4 

CM 

O 

X 

4 

4 

CD 

CM 

•—* 

O' 

X 

CM 

o 

4 

D 

O' 

4 

CM 

o 

X 

4 

o 

X 

iD 

ID 

-U 

O' 

4 

<D 

o 

X 

4 

4 

CM 

O 

X 

4 

CD 

4 

o 

4 

*'«- 

CD 

X 

4 

o 

4 

CM 

P*- 

D 

o 

X 

4 

O 

4 

CM 

CD 

O' 

d*' 

—H 

4 

CM 

X 

O 

i-U 

CM 

CM 

CD 

D 

4 

4 

4 

4 

4 

o 

• 

ID 

• 

• 

• 

CM 

• 

CM 

• 

PD 

• 

4 

• 

4 

• 

4 

• 

4 

• 

4 

• 

4 

• 

D- 

• 

X 

• 

X 

• 

O' 

• 

O' 

• 

• 

• 

•-i 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

CO  O  ^  OJ 

-«  i-H  ^  rsi  M  fsi 


Ti 


4*  r.  ^ 

<M  M  fsj 


«0  >  O  ^ 
fNj  <Ni  m  (o 


79 


THE  PCWIS  SPBCOKUM,  COTIREKCY,  AND  PHASE  REIATIONS  FOR  THE  NORTH-SOUTH  NACaiETIC  AND  EAST-WEST 
TELLURIC  FIELDS  AT  TBILISI  2  SEPTEMBER  1957  (LOW-PASS  FILTERED). 


In  addition  to  the  power  spectra,  Tables  9  through  l8  show  the  coherency 
(see  equation  56)  and  phase  angles  between  the  orthogonal  components  of  the 
telluric  and  magnetic  fields.  Figure  30  contains  a  few  plots  of  the  coherency 
between  the  components  of  the  magnetotelluric  fields*  l!bs  fbSM  ailfiM  teW 
an  erratic  behavior  and  hence  are  not  plotted. 

Comparison  with  other  Magnetotelluric  Power  Spectra 

Although  each  of  the  spectra  plotted  for  Lvov,  Ashkhabad,  and  Tbilisi 
cover  a  frequency  range  of  2.U  decades  of  frequency,  it  is  evident  that  this 
is  a  small  portion  of  the  total  magnetotelluric  power  spectrum.  In  order  to 
provide  a  more  comprehensive  viewpoint,  composite  spectra  are  plotted  in 
Figures  31  and  32  with  selections  from  data  computed  here  combined  with  spec¬ 
tra  obtained  from  other  authors. 

The  magnetotelluric  data  for  Ralston,  Canada,  were  recorded  on  31  July 
1959  and  were  supplied  by  C.  W.  Horton.  A  description  of  these  data  may  be 
found  in  papers  by  Horton  and  Hofftaan  [21,  22  j. 

The  magnetotelluric  data  for  Littleton,  Massachusetts,  were  recorded  \3(y 
Dr.  Thomas  Cantwell  on  5  December  1959*  These  spectra  are  given  in  his 
dissertation  [19]  as  Case  27-3* 

The  magnetotelluric  data  for  Austin,  Texas,  were  recorded  on  27  March  196I, 
by  Dr.  H.  W.  Smith  and  Mr.  F.  X.  Bostick  of  the  Electrical  Engineering  Research 
Laboratories,  Tiie  University  of  Texas.  Descriptions  of  the  recording  equip¬ 
ment  used  by  Smith  and  Bostick  are  found  in  the  reports  of  the  Electrical 
Bigineering  Research  Laboratories  [22,  23j- 

The  locations  of  these  various  stations  are  shewn  in  Figure  35*  Despite 
the  fact  that  the  data  were  recorded  over  a  period  of  2  l/2  years  and  extend 
over  a  longitude  range  of  170*,  the  data  appear  to  fit  a  uniform  trend.  It 
is  suggested  that  the  trends  established  in  Figures  and  32  are  a  first 
approximation  to  the  ambient  magnetotelluric  field  at  the  earth's  surface. 
Although  fluctuations  of  the  power  density  at  one  station  may  vary  10  db 
(Dr.  H.  W.  Smith,  personal  communication)  and  the  level  between  two  stations 
may  be  as  great  as  20  db  (Lvov  and  Tbilisi),  these  variations  are  small  com¬ 
pared  to  the  variation  of  dO  db  in  the  telluric  field  plotted  in  Figure  32* 
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FIGURE  30 

THE  COHERENCY  FUNCTION  Coh^^(f)  FOR  THE 
MAGNETOTELLURIC  FIELD  AT  TBILISI  1-2  SEPTEMBER  1957 
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FIGURE  31 

A  COMPOSITE  GRAPH  OF  THE  POWER  SPECTRA  OF  THE 
NORTH-SOUTH  MAGNETIC  FIELD  FOR  USSR,  CANADA, 
MASSACHUSETTS,  AND  TEXAS 
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FREQUENCY  IN  CYCLES  PER  KILOSECOND 
FIGURE  32 

COMPOSITE  GRAPH  OF  THE  POWER  SPECTRA  OF  THE 
EAST-WEST  TELLURIC  FIELD  FOR  USSR,  CANADA, 
MASSACHUSETTS,  AND  TEXAS 


The  corresponding  variation  In  the  x>ower  density  of  the  mapietlc  field 
plotted  In  Figure  31  Is  1^ 

The  following  empirical  formulas  represent  useful  approximations  to  the 
data  presented  In  Figures  31  32*  When  f  Is  expressed  in  <qreles  per 

kilosecond,  the  magnetic  field  can  he  expressed  by 


Pjj(f)  =  10^  1  +  25f^  +  f ^  +  1.3  X  10'*^  f^ 


-1-1 


(iMiMier/cps,(59) 


for  0.1  cpks  <  f  <  UOO  cpks. 


Superimposed  on  this  background  spectrum  Is  an  Intense  band  of  power  at 
80  cpks  In  the  Canadian  data.  Ihe  electric  field  can  be  expressed  by 
(again  f  In  cpks). 


©  /cps  ,  (60) 

for  0.03  cpks  <  f  <  400  cpks. 

The  experimental  data  for  P_(f)  scatter  more  strongly  than  the  data  for 
Pg(f).  This  Is  undoubtedly  due  to  the  Influence  of  the  earth's  structiire 
on  the  telluric  field. 


Pg(f)  =  10^ 


1  +  lOOf  + 


20f^ 


f5 
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CHAFTES  VII 


GEOPHYSICAL  INTERFRKIATIQN  OF  THE  PC^  SPECTRA 

Spectral  analysis  may  be  used  as  a  basis  for  geophysical  interpretation. 
Here  the  power  spectra  are  used  to  compute  an  apparent  resistivity  of  the 
earth  to  compute  the  frequency  response  of  the  earth  as  a  llneeur  filter,  and 
to  make  an  interesting  comparison  between  peaks  in  the  power  speetim  and 
modes  of  free  oscillation  of  the  earth. 

The  Apparent  Resistivity  of  the  Earth 

If  one  considers  the  energy  sources  of  the  magnetotelluric  field  to  be 
above  the  earth,  and  the  power  impinges  on  the  earth  as  a  plane  electromag¬ 
netic  wave,  most  of  the  power  is  reflected,  because  of  the  great  contrast 
in  the  electrical  properties  of  the  air  and  the  earth,  but  some  is  propagated 
into  the  earth.  To  a  first  approximation  the  tangential  electric  field  at 
the  earth's  surface  is  zero  and  the  tangential  magnetic  field  is  double  that 
of  the  i  -Ident  wave.  The  earth  is  not  a  perfect  conductor,  so  that  the 
tangential  electric  field  is  not  zero,  and  thus  is  sensitive  to  the  electri¬ 
cal  properties  of  the  earth.  The  reason  for  this  sensitivity  is  that  the 
incident  wave  and  the  waves  reflected  from  the  subsurface  strata  almost 
cancel  at  the  surface  [25].  The  extent  to  which  the  wave  incident  on  the 
eeurth  penetrates  the  earth  depends  upon  the  contrasting  electrical  proper¬ 
ties  of  a  layered  earth  and  the  power  spectrvun  of  the  incident  wave.  The 
depth  of  penetration  is  the  depth  at  which  the  magnitude  of  the  wave  is 
equal  to  l/e  of  its  magnitude  at  the  surface.  This  depth  is  called  the  skin 
depth  and  is  a  function  of  frequency  [26]. 

Assvme  that  the  air  and  earth  are  each  semi -infinite  electrically  homoge¬ 
neous,  isotropic  media  such  that  the  air-earth  Interface  is  a  plane  located  at 
z  =  0  in  the  right-handed  rectangular  cartesian  coordinate  system  shown  in 
Figure  54. 

The  electric  and  magnetic  fields  at  any  point  in  a  homogeneous,  laotroplc, 
euid  charge  free  medliuD  satisfy  Maxwell's  equations  in  the  form  [26] 


FIGURE  34 

THE  AIR- EARTH  INTERFACE 
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(6la) 

(6lb) 

(6lc) 

(6ld) 


7  X  E  =  5:^  ' 

Vxi?=-^+6^ 

p  ot  ^ 

V  •  ^  =  0  , 

and  7  *  if  =  0  , 

where  ^  is  the  electric  field  strength  in  volts  per  meter, 

^  is  the  magnetic  field  strength  in  amperes  per  meter, 

H  is  the  magnetic  permeability  in  henry  per  meter, 
e  is  electric  permittivity  in  farad  per  meter, 
p(=  is  the  resistivity  in  ohm-meters,  and 

a  is  the  conductivity  in  mho  per  meter. 

The  usual  conditions  for  neglecting  the  displacement  current  [26],  i.e. 


r^^  »i 


(62) 


are  amply  met  here  so  that  equation  (611))  may  be  approximated  by 

7  X  ^  i  ^  . 

p 


(63) 


If  it  be  assumed  that  linearly  polarized  plane  electromagnetic  waves 
are  propagated  in  the  positive  z  direction  such  that  the  x-axis  is  in  the 
direction  of  the  electric  field  and  the  y-axis  is  in  the  direction  of  the 
magnetic  field,  then  there  is  no  variation  of  the  fields  in  the  x  or  y 
directions.  For  the  above  conditions,  equations  (6la)  and  (63)  reduce  to 


Se 

x 


= 


6h 


(64) 


and 


6h 

=  -0  E 

OZ  X 


(65) 
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where  the  subscripts  x  and  y  Indicate  the  polarization.  Substitution  of 
equation  (65)  into  equation  (64)  or  equation  (64)  into  (65)  yield 


Sh 

=  5t^ 


oz 

since  the  fields  will  be  resolved  in  terms  of  harmonic  components  it  is 
natural  to  select  harmonic  solutions  of  the  form 

E(z,(a,t)  =  E^(z,a))  e*^^^  ,  (68) 


Hy(z,cu,t)  =  Hy(z,aj)  e*^  .  ( 

Substituting  equations  (68)  and  (69)  into  (64)  throu^  (67)  results  in 

8e  (zjO)) 

- ^ -  =  -Jqj  Hy(z,a3)  ,  ( 

8h  (z,a>) 


S  H  (z,co) 

- *-5 -  =  Jcuou  H  (z,cu)  , 


8  E^(z,a3) 


=  juxjp  E^(z,co) 
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Solutions  of  equations  (72)  and  (75)  can  be  written  in  the  form 


Hy(z,a))  =  H^y(a))  e"J^  ,  (74) 

and 

E^(z,a))  =  E^^(a))  e'’^  ,  (75) 

where 

=  -jooii  ,  (76) 

or 


k  = 


(1-j) 


(77) 


Introducing  the  solutions  (74)  and  (75)  into  the  component  equation  (65) 
results  in 


E  (c))  .. 

ox _  _  jk 

H  (cu)  ""  a 
oy' 


or 


E  (o)) 


ox 


iTTSi 

oy' 


Using  the  identity 

l+J  .  eJ"/**  .  ■/?  , 

and  solving  for  the  resistivity^  eqviation  (78)  becomes 


P  = 


2rtuf 


E  (f) 

OX^ 


H  (fj 


oy 


where  f  = 


(78) 


(79) 

(80) 
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Assuming  that  the  magnetotellurlc  field  Is  due  to  a  linearly  polarized 
plane  electromagnetic  wave  propagating  from  the  air  Into  the  earth  (see 
Figure  3^)  It  is  easy  to  see  how  the  resistivity  of  a  homogeneous,  isotropic 
earth  could  be  computed  from  the  magnetotellurlc  field  components  using 
equation  (80),  provided  one  could  separate  the  frequency  components  in  the 
experimental  observations. 

Cagniard's  paper  [24]  on  the  magnetotellurlc  field  produced  a  renewal 
of  interest  in  the  use  of  the  magnetotellurlc  field  components  for  the 
interpretation  of  the  electrical  8ti*uctiu*e  of  the  earth.  Cagnlard  assumes 
stratified  earth  models  in  which  the  earth  consists  of  two  or  more  horizontal 
layers  of  contrasting  resistivity.  Equation  (8o),  which  gives  the  resistivity 
for  an  unlayered  earth,  may  also  be  used  for  a  layered  earth,  in  which  case 
the  computed  quantity,  p,  is  referred  to  as  the  apparent  resistivity,  p  . 

St 

Cagniard*8  proposal  involves  the  plotting  of  an  apparent  resistivity  versus 
the  square  root  of  the  period  of  the  wave,  hut  he  does  not  show  how  one  should 
resolve  the  experimental  observations  which  resemble  a  stochastic  time  series 
into  frequency  components. 

Cantwell  considered  the  magnetotelluric  field  components  as  samples  of 
stationary  stochastic  processes  and  applied  the  numerical  methods  of  Blackman 
and  Tukey  to  the  data.  Cantwell  showed  that  Cagniard*s  basic  formxila  for  the 
apparent  resistivity  may  be  expressed  in  terms  of  power  spectra  as 

Pp(f) 

».  ■  ^  '  <81) 

n 

where  p  =  apparent  resistivity,  ohm  meters, 
a 

T  =  period  in  seconds, 

Pg(f)  =  power  spectral  density  of  the  telluric  field  component  in  (milli- 
volts/kilometer)  per  cycle  per  second, 

Pjj(f)  =  power  spectral  density  of  the  orthogonal  component  of  the  horizon- 
tal  magnetic  field  in  (gammas)  per  cycle  per  second, 

0.2  =  dimensional  constant  chosen  to  make  the  equation  correct. 
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Cagnlard  provides  a  template  for  interpreting  the  apparent  resistivity 
in  terms  of  electrically  conducting  layers  when  the  data  are  plotted  as  a 

1/2 

gjraph  of  p  versus  T  '  on  log-log  paper.  However,  the  present  data  do  not 
a 

show  evidence  of  an  electrical  interface  so  they  are  plotted  versus  frequency. 
The  vertical  scale  is  compressed  to  reduce  the  scatter  of  the  points  and  to 
facilitate  the  drawing  of  a  curve  throu^  the  points.  Figure  35  contains 
plots  of  apparent  resistivity  hased  on  the  Bast -West  telluric  field  and  the 
North-South  magnetic  field  for  1  and  2  September.  Figure  36  contains  a 
similar  plot  based  on  the  other  pair  of  components  of  the  magnetotelluric 
field.  The  computed  values  of  apparent  resistivity  are  given  in  Tables  19-22. 
The  latter  data  are  restricted  to  the  outputs  of  the  high-pass  filter  because 
of  the  uncertainties  in  the  low -passed  data. 

The  four  curves  plotted  in  Figures  35  and  36  are  mutually  consistent 
and  they  indicate  a  resistivity  of  6  ohm  meters  over  the  frequency  range  0.1 
to  T  cpks.  The  apparent  resistivity  is  essentially  constant  over  this  fre¬ 
quency  range  and  there  is  no  evidence  of  stratification  of  the  earth's 
conductivity.  The  difference  between  the  apparent  resistivities  obtained 
from  the  two  components  of  the  telluric  field  are  not  sufficiently  different 
to  suggest  any  horizontal  anisotropy  in  the  earth's  resistivity  at  the  rele¬ 
vant  depths - 

This  determination  of  the  apparent  resistivity  of  the  earth  is  not  incon¬ 
sistent  with  the  results  of  previous  theoretical  and  experimental  studies. 
Cantwell  interprets  the  resistivity  profile  in  Massachusetts  as  a  layer  70 
kilometers  thick  with  a  resistivity  of  8,000  ohm  meters  underlain  by  a  zone 
of  much  lower  resistivity.  The  frequency  range  covered  by  the  present  study 
is  not  so  low  that  the  layer  determined  by  Cantwell  has  no  effect  on  the 
measurements.  However,  the  results  of  Figures  35  and  36  suggest  that  the 
resistivity  of  the  material  below  the  layer  determined  by  Cantwell  is  6  ohm 
meters.  The  present  data  do  not  establish  a  thickness  for  either  of  the 
layers  nor  of  the  resistivity  contrast.  However,  since  a  change  of  frequency 
of  a  factor  of  100  is  equivalent  to  a  change  of  depth  of  a  factor  of  10  for 
resistivity  determinations  [25] >  one  feels  safe  in  concluding  that  the  layer 
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APPARENT  RESISTIVITY  APPARENT  RESISTIVITY 

IN  OHM-METERS  _  IN  OHM-METERS 


THE 


FIGURE  35 

APPARENT  RESISTIVITY  VERSUS  FREQUENCY  DETERMINED 
FROM  EAST-WEST  COMPONENT  OF  TELLURIC  FIELD 
AT  TBILISI  I  AND  2  SEPTEMBER  1957 
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FREQUENCY  IN  CYCLES  PER  KILOSECOND 
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FREQUENCY  IN  CYCLES  PER  KILOSECOND 


FIGURE  36 

THE  APPARENT  RESISTIVITY  VERSUS  FREQUENCY  DETERMINED 
FROM  NORTH-SOUTH  COMPONENT  OF  TELLURIC  FIELD 
AT  TBILISI  I  AND  2  SEPTEMBER  1957 
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APPARENT  RESISTIVITY  IN  OHM  METER 


FROUENCY 

PERIOD 

EAST  MAGNETIC 

NORTH  MAGNETIC 

CPS 

SEC 

NORTH  TELLURIC 

EAST  TELLURIC 

1 

•OOOE+00 

.00 

.OOOOOOE+00 

.OOOOOOE+00 

2 

.231E-03 

4320.00 

.369279E+01 

.281286E-«-02 

2 

.463E-03 

2160.00 

.113957E+02 

.194253E+02 

4 

.694E-03 

1440.00 

.840330E+01 

.129704E+02 

c 

.926E-03 

loco. 00 

.521944E+01 

. 144856E+02 

6 

.116E-02 

864.00 

.436688E+01 

.129919E+02 

7 

.139E-02 

720.00 

.531387E+C1 

.600730E+C1 

8 

.162E-02 

617.14 

.427588E+01 

.463612E+01 

9 

.185E-02 

540.00 

.519203E+01 

.465192E+01 

10 

.208E-02 

480.00 

.523871E+01 

.343634E+01 

11 

.231E-02 

432.00 

.252099E+01 

.265340E+01 

12 

.255E-02 

392.73 

.220273E+01 

.234180E+01 

13 

.278E-02 

360.00 

.319496E+01 

.298110E+01 

14 

.301E-'02 

332.31 

.448339E+01 

.391637E+01 

15 

.324E-02 

308.57 

.575062E+01 

.632859E+01 

16 

.347E-02 

288.00 

.829431E+01 

.548903E+01 

17 

.370E-02 

270.00 

.135819E+02 

.461995E+01 

18 

.394E-02 

254.12 

.163477E+02 

.820012E+01 

19 

.417E-02 

240.00 

.129818E+02 

.948598E+01 

20 

.440E-02 

227.37 

.108044E+02 

.337316E+01 

21 

.463E-02 

216.00 

.829462E+01 

.242292E+01 

2? 

.486E-02 

205.71 

.608538E+C1 

.215063E+C1 

2- 

.509E-02 

196.36 

.834416E+C1 

. 190608E+01 

24 

.532E-02 

187.83 

.110422E+02 

.273064E+C1 

25 

.556E-02 

180.00 

.802813E+01 

.533159E+01 

26 

.579E-02 

172.80 

.717575E+C1 

. 145645E+C2 

27 

.602E-02 

166.15 

.811710E+01 

.227996E+02 

28 

.625E-02 

160.00 

.927438E+01 

.103059E+02 

29 

.648E-02 

154.29 

.101736E+02 

.316596E+01 

30 

•671E-02 

148.97 

.961852E+01 

.494256E+01 

31 

.694E-02 

144.00 

.821618E+01 

.621548E+01 

TABLE  19 

APPARENT  RESISnVITy  AT  TBILISI  1  SEPTEMBER  1957  (HIGH-PASS  FILTERED) . 
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APPARENT  RESISTIVITY  IN  OHM  METER 


FROUENCY  PERIOD 

CPS  SEC 


EAST  MAGNETIC  NORTH  MAGNETIC 
NORTH  TELLURIC  EAST  TELLURIC 


1 

•OOOE+00 

.00 

2 

.579E-04 

17280.00 

3 

.116E-03 

8640.00 

4 

.174E-03 

5760.00 

5 

.231E-03 

!^320.00 

6 

.289E-03 

3456.00 

7 

.347E-03 

2880.00 

8 

.405E-03 

2468.57 

9 

•463E-03 

2160.00 

10 

.521E-03 

1920.00 

11 

.579E-03 

1728.00 

12 

.637E-03 

1570.91 

13 

•694E-03 

1440.00 

14 

.752E-03 

1329.23 

15 

.810E-03 

1234.29 

16 

.868E-03 

1152.00 

17 

.926E-03 

1080.00 

18 

.984t-03 

I016r47 

19 

.104E-02 

960.00 

20 

.llOE-02 

909.47 

21 

.116E-02 

864.00 

22 

.122E-02 

822.86 

23 

.127E-02 

785.46 

24 

.133fc-02 

751.30 

25 

.139E-02 

720.00 

26 

.146E-02 

691.20 

27 

.150E-02 

664.62 

28 

.156E-02 

640.00 

29 

.162E-02 

617. 14 

30 

.168E-02 

595.86 

31 

.1 74E-02 

5  76.00 

•OOOOOOE+00 
• 500494E+00 
-.123966E+02 
-.327064E+02 
-.185662E+02 
.449952E+02 
•822454E+02 
.592236E+02 
-.:>35531E+02 
-.811736E+02 
-•  177620E+02 
.129386E+02 
-.386669E+02 
-.501727E+02 
-.701871E+01 
.3161O0E+O2 
.937247E+02 
.162591E+03 
-.281032E+02 
-.588788E+02 
.145  700E+02 
.238263E+02 
-.234684E+02 
-.656Vb2E+02 
-.46b601E+02 
-.33  7360E  +  02 
-.582475E+02 
-.120231E+02 
.264310E+02 
. 149784E+02 
.  7808b5E-<-0i 


.OOOOOOE+00 
.460016E+01 
.912174E+01 
.794878E+02 
.970748t+0l 
.318607E+02 
.132243E+02 
.859154E+01 
.109381E+02 
.135725E+02 
.103109E+02 
.665483E+01 
.344252E+01 
.294153E+01 
.261474E4-01 
.241b98E+01 
. 104  742E+02 
-.108015E+02 
“.168648E+02 
•164810E+03 
• 710142E+01 
•471713h+01 
.675l95t+01 
.6b043/E+01 
.8^8513£+01 
.2306V3E+02 
•149382E+02 
.6062 1 3  E+0 1 
•2228616+01 
•810969E+00 
•602386E+00 


TABLE  20 

APPARENT  RESISTIVITY  AT  TBILISI  1  SEPTEMBiR  1957  (LOW-PASS  FILTERED). 
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APPARENT  RESISTIVITY  IN  OHM  METER 


FROUENCY 

PERIOD 

EAST  MAGNETIC 

NORTH  MAGNETIC 

CPS 

SEC 

NORTH  TELLURIC 

EAST  TELLURIC 

1 

•OOOE+00 

•  00 

.OOOOOOE+00 

•OOOOOOE+OO 

2 

.231E-03 

4320*00 

•848561E^00 

.705163E+01 

3 

.463E-03 

2160.00 

.112731E+01 

.720980E+01 

4 

•694E-03 

1440.00 

•  208879E-»-01 

•  805363E-I-01 

5 

.926E-03 

1080.00 

•256123E+01 

.719720E+01 

6 

.116E-02 

864.00 

•  222840E-I-01 

.596853E+01 

7 

•139E-02 

720.00 

.240491E+01 

.546525E+01 

8 

.162E-02 

617.14 

.243540E+01 

.613543E+01 

9 

.185E-02 

540.00 

•  201286E-^01 

.511269E+01 

10 

.208E-02 

480.00 

.186138E+01 

.414423E+01 

11 

.231E-02 

432.00 

.200885E+01 

.470676E+01 

12 

.255E-02 

392.73 

.222912E+01 

.534337E+01 

13 

.278E-02 

360.00 

.219383E+01 

•586071E+01 

14 

.301E-02 

332.31 

.219664E+01 

.664984E+01 

15 

.324E-02 

308.57 

.260272E+01 

.635374E+01 

16 

.347E-02 

288.00 

•581040E+01 

•677598E+01 

17 

.370E-02 

270.00 

.996367E+01 

•587172E+01 

18 

.394E-02 

254.12 

.583734E+01 

•  520250E-*'0l 

IS 

.417E-02 

240.00 

.335621E+01 

•641130E+01 

20 

.440E-02 

227.37 

•399993E+01 

•877473E+01 

21 

.463E-02 

216.00 

.443690E+01 

.707639E+01 

22 

.486E-02 

205.71 

.446603E+01 

•762637E+01 

23 

.509E-02 

196.36 

,444960E+01 

•862782E+01 

24 

.532E-02 

187.83 

.539938E+01 

•867567E+01 

25 

.556E-02 

183.00 

.477207E+01 

.558615E+01 

26 

.579E-02 

1"2.80 

.349082E+01 

.432838E+01 

27 

.602E-02 

166.15 

•374401E+01 

.563683E+01 

28 

.625E-02 

160.00 

.292866E+01 

.586957E+01 

29 

.648E-02 

154.29 

.209665E+01 

.433147E+01 

30 

.671E-02 

148.97 

.198033E+01 

•  568006E-I-01 

31 

.694E-02 

144.00 

.202499E+01 

•  831898E-»-0l 

TABLE  21 

APPARENT  RESISnViry  AT  TBIUSI  2  SEPTEMBER  1957  (HIGH-PASS  FILTERED) 
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APPARENT  RESISTIVITY  IN  OHM  METER 


FRQUENCY 

PERIOD 

EAST  MAGNETIC 

NORTH  MAGNETIC 

CPS 

SEC 

NORTH  TELLURIC 

EAST  TELLURIC 

1 

•OOOE+00 

•  00 

•  OOOOOOE-^OO 

•OOOOOOE+00 

2 

•579E-04 

17280.00 

-•450131E+01 

•753286E+01 

3 

•116E-03 

8640*00 

.515451E+02 

.622565E+01 

4 

.174E-03 

5760.00 

.108739E+03 

•134203E+02 

5 

.231E-03 

4320.00 

-.399511E+03 

.631158E+01 

6 

•289E-03 

3456.00 

-•J42292E+03 

.633122E+01 

7 

.347E-03 

2880.00 

.290387E+03 

.696730E+01 

8 

•405E-03 

2468.57 

.152132E+03 

•102615E+02 

9 

.463E-03 

2160.00 

-•160802E+03 

.105180E+02 

10 

.521E-03 

1920.00 

•115600E+01 

.154988E+02 

11 

.579E-03 

1728.00 

.604839E+03 

.652822E+01 

12 

•637F-03 

1570.91 

.503195E+03 

.475916E+01 

13 

.694E-03 

1440.00 

.926298E+01 

.468358E+01 

14 

•752E-03 

1329.23 

-.176087E+03 

.422466E+01 

15 

.810E-03 

1234.29 

-.157443E+03 

.318379E+01 

16 

.868E-03 

1152.00 

-.145195E+03 

.300343E+01 

17 

•926E~03 

1080.00 

-•.550225E+02 

.569165E+01 

18 

.984E-03 

1016.47 

-.947606E+03 

.951582E+01 

19 

.104E-02 

960.00 

-•334956E+04 

.430453E+01 

20 

.110E~02 

909.47 

-.107317F+04 

.393863E+01 

21 

.116E~02 

864.00 

.691856E+01 

.440579E+01 

22 

.122E~02 

822.86 

-.752885E+02 

.731142E+01 

23 

.127E~02 

785.45 

-.232905E+03 

.131430E+02 

24 

.133E~02 

751.30 

-.482368E+03 

.536235E+01 

25 

.139E-02 

720.00 

-.148510E+C3 

.529453E+01 

26 

.145E-02 

691.20 

.376149E+03 

.967943E+01 

27 

•150E~02 

664.62 

.520689E+03 

.197032E+02 

28 

.156E-02 

640.00 

-.591264E+02 

.372136E+62 

29 

•162E-02 

617.14 

-.131727E+03 

.634290E+01 

30 

•168E-02 

595.86 

-.270762E+91 

.150332E+01 

31 

.174E-02 

576.00 

•  124518E-t-02 

.110285E-F01 

TABLE  22 

APPARENT  RESianVITT  AT  TBILISI  2  SEPTEMBER  1957  (LCW-PASS  Pn/PIRED) 
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whose  resistivity  is  6  ohm  meters  extends  from  about  70  to  700  kilometers 
below  the  surface. 

The  survey  paper  by  Tozer  [l]  shows  plots  of  resistivity  versus  depth 
computed  theoretically  by  several  authors.  These  plots  indicate  that  a 
resistivity  of  6  ohm  meters  is  typical  of  a  depth  of  500  kilometers. 


Freauenca 


3nse  of  the  Barth 


Use  of  Random  Functions  to  Determine  the 


Frequency  Response  of  a  Linear  Filter 


Consider  a  linear  filter  such  as  shown  in  Figure  57-  Let  the  input, 
u(t),  be  a  member  function  of  a  stationary  random  process.  The  output,  v(t), 
will  also  be  a  stationary  random  process  [9,  I8].  Recalling  equation  (l), 
the  output  may  be  written  in  terms  of  the  input  throu^  the  relation 


p  00 

'■(t)  =  /  u(t-a)  W(a)  da  , 


where  W  is  the  impulsive  response  of  the  filter.  Shifting  the  time  scale  by 
an  amount  t  results  in 

p  00 

v(t+T)  =  /  u(t+T-a)  W(a)  do  .  (85) 

Multiplying  (83)  by  u(t)  and  taking  u(t)  inside  the  integral  results  in 


u(t)  v(t+T)  =  /  u(t)  u(t+T-a)  W(a)  da  .  (84) 

Integrating  with  respect  to  t,  multiplying  by  1/t,  and  taking  the  limit  as  t 
increases  without  limit  results  in 

^ ^  LL jL  L- 


...(85) 


The  order  of  integration  on  the  right  hand  side  of  equation  (85)  may  be 


inverted  so  that 


100 


llm 


h  f  ^  u(t)  u(t+T)  dT  =  r  w(a)  f  u(t)  u(t+T-a)  dTda 

T  J.T/2  ^0"  ^-T/2 


T/2 


...(86) 

Using  the  definitions  of  auto-  and  cross -correlations  (86)  may  he  written  as 

C^v^t)  =  /^"w(a)  C^Jt-c)  dc  .  (87) 

IJaking  the  Fourier  transform  of  (87)  and  using  the  definition  of  Y(f)  resvilts  . 
in 


P  (f)  =  Y(f)  P  (f) 
uv'  '  uu' 


(88) 


or 


P  (f) 
uu' 


The  phase  delay  is  given  by 


0(f)  =  arc  tan 


p°'*^  (f) 

^uv  ' 


pCven 

uv 


(f) 


(89) 


(90) 


in  accordance  with  the  previous  discussion  of  cross-power.  Equation  (89) 
shows  how  the  frequency  response  of  a  linear  filter  may  be  determined  from  the 
ratio  of  the  cross-power  (between  the  input  and  output)  to  the  auto  power  of 


the  input. 


Frequency  B^sponse  of  the  Earth  as  a  Linear  Filter 


Assuming  that  the  earth  behaves  as  a  linear  filter  with  the  variation 
of  a  horizontal  component  of  the  earth's  magnetic  field  as  the  random  input 
or  driving  function  and  the  corresponding  component  of  the  telluric  field  as 
the  output  or  resulting  function,  the  frequency  response  of  the  earth  can 
be  computed  using  equation  (89). 
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Specifically,  the  frequency  response  of  the  earth  is  computed  from 


where 


^earth^^^  "  Pjjjj(f) 


(91) 


is  the  magnitude  of  the  cross -power  between  the  magnetic  and 
electric  field  components, 


P  (f)  is  the  (auto-)  power  spectrum  of  the  magnetic  field  component. 

lUl 

Results  using  the  above  relation  are  tabulated  in  Tables  23,  24,  25,  and  26, 
and  plotted  in  Figures  58  and  39*  These  plots  show  that  the  frequency 
response  is  proportional  to  the  square  root  of  frequency,  a  characteristic  of 
a  uniform  conductor. 

It  is  felt  that  (91)  is  a  more  significant  way  of  looking  at  the 
electrical  properties  of  the  earth  since  it  includes  the  coherency  between 
the  components  whereas  Cagniard's  formula  as  used  by  Cantwell  (81)  does  not. 


A  Comparison  Between  Observed  Free  Earth  Periods 
and  Peaks  in  the  Geomagnetic  Power  Spectrum 

A  recent  paper  by  Ness,  Harrison,  and  Slichter  [25]  reports  values  of 
the  periods  of  the  free  oscillation  of  the  earth  as  observed  with  a  gravimeter 
during  the  Chilean  earthquake  of  22  May  i960.  These  reported  values  of  the 
periods  of  the  spherical  modes  are  of  interest  here  since  the  frequency  range 
of  the  power  spectra  computed  from  low-pass  filtered  data  recorded  at  Tbilisi 
is  the  same  as  the  frequency  range  of  the  oscillations  of  the  earth. 

One  would  expect  that  the  earth  (a  conductor)  oscillating  in  its  magnetic 
field  would  experience  induced  alternating  currents  with  frequencies  corre¬ 
sponding  to  the  frequency  of  oscillation.  With  this  in  mind,  the  data  for 
1  and  2  September  were  subjected  to  power  spectral  analysis  using  60  lags 
instead  of  30  lags  as  were  used  to  compute  the  results  shown  in  Chapter  IV. 

The  results  of  this  computation  are  plotted  in  Figures  40  and  4l.  To  permit 
careful  examination  of  all  peaks  in  the  power  spectra,  the  graphs  were  plotted 
on  semi -logarithmic  paper  with  frequency  on  the  linear  scale.  Table  27  shows 
a  comparison  between  the  observed  free  oscillation  modes  and  the  peaks  in  the 
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FREQUENCY  RESPONSE  OF  THE  EARTH 


FROUENCY 

PET lOD 

EAST  MAGNETIC 

NORTH  MAGNETIC 

CPS 

SEC 

NORTH  TELLURIC 

EAST  TELLURIC 

1 

•OOOE+00 

•  00 

.OOOOOOE+00 

.OOOOOOE+00 

2 

.231E-03 

4320.00 

.101969E-01 

.853222E-01 

2 

•463E-03 

2160*00 

.493573E-01 

.149597E+00 

A 

.694E~03 

1440.00 

.669767E-01 

.123002E+00 

2 

.926E-03 

1080.00 

.302052E-01 

.104105E+00 

6 

.116E-02 

864.00 

.631482E-01 

.140551E+00 

7 

•139E-02 

720.00 

.346207E-01 

.121555E+00 

8 

•162E-02 

617.14 

.920681E-01 

.733017E-01 

9 

.185E-02 

540.00 

.127906E+00 

.121134E-01 

10 

.208E-02 

480.00 

.111542E+00 

.162631E-01 

11 

.231E-02 

432.00 

.346192E-01 

.240620E-01 

12 

.255E-02 

392.73 

.325292E-01 

.352688E-02 

13 

.278E-02 

360.00 

.600825E-01 

.245335E-01 

14 

.301E-02 

332.31 

.119530E+00 

.328242E-01 

15 

.324E-02 

308.57 

.120579E+00 

.642141E-01 

16 

.347E-02 

288.00 

.492887E-01 

.914414E-01 

17 

.370E-02 

270.00 

.187579E+00 

.583060E-01 

18 

.394E-02 

254.12 

.229296E+00 

.897950E-01 

19 

.417E-02 

240.00 

.185926E+00 

.138230E+00 

20 

•440E-02 

227.37 

.100954E+00 

.501925E-01 

21 

.463E-02 

216.00 

.114234E+00 

.104349E+00 

22 

.486E-02 

205.71 

.153928E+00 

.616915E-01 

23 

.509E-02 

196.36 

.177725E+00 

.326258E-01 

24 

.532E-02 

187.83 

.174232E+00 

.118238E-01 

25 

.556E-02 

180.00 

.695956E-01 

.449745E-01 

26 

•579E-02 

172.80 

.668888E-01 

.104086E+00 

27 

.602E-02 

166.15 

.142239E+00 

.113472E+00 

26 

.625E-02 

160.00 

.120917E+00 

.519344E-01 

29 

.648E-02 

154.29 

.126937E+00 

.668369E-01 

30 

.671E-02 

148.97 

.110482E+00 

.873261E-01 

31 

.694E-02 

144.00 

.126772E+00 

.545266E-01 

TABLE  25 

FREQUENOY  RESPONSE  OF  THE  EARTH  AT  TBILISI  1  SEPTEMBER  1957 
(HKZ-FASS  filtered). 
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0  <30  T'  o  ^  ^  Xf\  o  I <33  -  A  o  ►-<  rg  rrv  li^  O  ^-  tO  CA  c>  •«! 

.-•  »H  f-i  f-H  .H  *  4  r>j  I  <4  r4  r  4  rj  r j  c  g  c  4  ( >i  t  j  m 


fRlQUL.ti'Zi  RLGPONoE  OF  THE  EARTH 


frquency 

F  wA i ww 

EAST  MAGNETIC 

NORTH  MAGNETIC 

CPS 

SlC 

NORTH  TELLURIC 

EAST  TELLURIC 

•QOOE+00 

•  00 

•OOOOOOE+00 

.0000002+00 

.579E-04 

17280.00 

. 6118772-01 

.19939*fE-0l 

.116L-03 

8640.00 

.  146646E  +  00 

.4723942-01 

.i74E-03 

5760.00 

.526207E+00 

.1745142+00 

.231E  -03 

4320.00 

.2842712+00 

•  7 125512-01 

.2SQL-03 

3456.00 

.552774E+C0 

.1139322+00 

•  34  r  tl-03 

2880.00 

.6925952^00 

.0993332-01 

.405E-03 

2468 .57 

.7924412+00 

.9535342-01 

.•+b3L-03 

2160.00 

.984499L+00 

.1352642+00 

.321t:-03 

1920.00 

.1201302+01 

.1337822+00 

•  3  7  9 1_  0  3 

1  (28.00 

•  680  756L  +  00 

.9h47332-0; 

•637L '03 

1570.51 

.2846792+00 

.6163332-01 

•  t>34i.-03 

i.  440 .00 

•236C46L+00 

.4  734942  -Cl 

.752F-03 

1329.23 

.1931892+00 

.3078852-01 

•  &  i.  0  C.  '  yj  Jt 

1234.29 

•  2362  78E  +  00 

.1389902-01 

.368E-03 

1152.00 

.429752E+00 

.2275542-01 

•  72&C.  03 

lOSO.Ov; 

•  103  (’162  +  01 

.9584682-01 

.364E  03 

1016.47 

•35840 1 2+01 

.2943142-01 

•  1 04 1.  0 

960.00 

•  2377  55L  +  Ga 

•676752L-01 

.llCE-02 

909.47 

.341329E+00 

•  143 1 06E  +  0 1 

.116L  02 

8o4.00 

•37 9 o 13u+00 

.6961562-01 

•  122E-02 

62  2 • 86 

.838526E+00 

.1633152-01 

•  X  2  7 1_  02 

1  0  3  •  4  5 

•5759552+00 

•  A  72  7482  ♦'00 

.1332-02 

7:  1.30 

.5852122-01 

•  1  76i34c.  +  00 

•1392- 02 

('20.00 

•  79  5 1 062  +  00 

.5966012-01 

.1452-02 

691.20 

.li-r933E  +  0i 

.1316062+00 

.1502  -02 

o6  4 .62 

•2659i3L*0i 

.5H205i2-0i 

.13GE-02 

64  0 • GO 

.2539022<0i 

.1513272+00 

.1622  0^. 

617.14 

•13692 5Li0i 

•  1202612+00 

.168E-02 

595.36 

.^124162+00 

.5207092-01 

.174(.  02 

576.00 

.156070L+00 

.20  79172-01 

TABLE  2k 

FREQUENCY  RESPONSE  OF  THE  EARTH  AT  TBILISI  1  SEPTEMBER  19^7 
(LON-PASS  FILTERED) 
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FREQUENCY  RESPONSE  OF  THE  EARTH 


FRQUENCY 

CPS 


PERIOD 

SEC 


EAST  MAGNETIC  NORTH  MAGNETIC 
NORTH  TELLURIC  EAST  TELLURIC 


1 

•OOOE+OO 

•  00 

2 

.231E-03 

4320*00 

3 

.463E-03 

2160*00 

4 

.694E-03 

1440*00 

5 

.926E-03 

1080*00 

6 

.116E-02 

864*00 

7 

.139E-02 

720*00 

8 

.162E-02 

617*14 

9 

.185E-02 

540*00 

1C 

.208E-02 

480*00 

11 

.231E-02 

432.00 

12 

.255E-02 

392*73 

13 

.278E-02 

360.00 

14 

.301E-02 

332*31 

15 

.324E-02 

308*57 

16 

.347E-02 

283*00 

17 

.370E-02 

270*00 

18 

.394E-02 

254*12 

19 

.417E-02 

240*00 

20 

.440E-02 

227*37 

21 

.463E-02 

216*00 

22 

.486E-02 

205*71 

23 

.509E-02 

196*36 

24 

.532E-02 

187.83 

25 

.556E-02 

180.00 

26 

.579E-02 

172.80 

27 

.602E-02 

166.1 5 

28 

.625E-02 

160.00 

29 

.648E-02 

154*29 

30 

.671E-02 

148*97 

31 

.694E-02 

144*00 

•OOOOOOE+OO 

.976764E-02 

.910998E-02 

.962849E-02 

.155341E-01 

.257682E-01 

.505343E-01 

.636224E-01 

.508070E-01 

.697188E-01 

.812226E-01 

.552767E-01 

.638727E-01 

.530786E-01 

.467134E-01 

.559962E-01 

.141901E+00 

.105545E+00 

.416303E-01 

.162070E-01 

.461712E-01 

.981053E-01 

.225629E-01 

.828101E-01 

.739639E-01 

.787050E-01 

.934917E-01 

.782442E-01 

.640537E-01 

.572892E-01 

.550426E-01 


•OOOOOOE+OO 

.650262E-01 

.948715E-01 

•928413E-01 

•101237E+00 

•125262E+00 

.121177E+00 

.877263E-01 

•539094E-01 

.980794E-01 

.880815E-01 

.995669E-01 

.148124E+00 

•172489E+00 

.131634E+00 

• 105226E+00 

•113286E+00 

.162493E+00 

.142467E+00 

•887039E-01 

•179781E+00 

.141287E+00 

.929740E-01 

.157651E+00 

.115606E+00 

.857313E-01 

• 162163E+00 

.178053E+00 

.122874E+00 

.119543E+00 

•141179E+00 


TABLE  25 

FREQUENCY  RESPONSE  OF  THE  EARTH  AT  TBILISI  2  SEPTEMBER  1957 
(HIC2I-PASS  FILTERED). 
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FREQUENCY  RESPONSE  OF  THE  EARTH 


FROUENCY 

PERIOD 

EASF  MAGNETIC 

NORTH  MAGNETIC 

CPS 

i£C 

NORTH  TEuLURIC 

EAST  TELLURIC 

1 

•OOOE+OO 

•  00 

.OOOOOOE+00 

•OOOOOOE+00 

2 

•Sy9E-04 

1  T2«0.00 

•  73318  7E-01 

.337248E-01 

3 

•116E-03 

8640.00 

•  6  76029E  +  00 

.603887E-01 

4 

•174E-03 

3760.00 

•347113E+01 

.938761E-01 

6 

.231E-03 

4320.00 

.41U314E+01 

.812649E-01 

6 

•289E-03 

3436.00 

.329376E+01 

•904930E-01 

7 

•347E-03 

2880.0J 

•  19488  7E  +  01 

•104736E+00 

8 

.40SE-03 

2468.3  / 

.86632  7E+01 

. 1 14438E+00 

9 

.463E-03 

2160.00 

.8336  /lE  +  01 

•104173E+00 

10 

.52iE-03 

1920.00 

• 1 1 3806E+02 

.110836E+00 

11 

•  S  /yh-03 

1  /28.00 

•  1 1  2466E-*-02 

•9d0232E-0l 

12 

.63  TE-03 

1370.91 

•  8  7  7U49E  +  01 

• 108320E+00 

13 

.694E-03 

1440.00 

.389231E+01 

.92360bE-0l 

14 

• /S2E-03 

1329.23 

•233444E+01 

.697039E-01 

l^> 

•BlOE-03 

1234.29 

•  18  u260E  +  0 1 

•  70884  7E-01 

16 

•66«E-03 

1132.00 

•241089E+01 

• 703739E-01 

1  / 

•v2bt-03 

1080.00 

•8 7063UE+Oi 

•932386t-01 

lo 

•9b4E-03 

1016.47 

• i61v06E+02 

.228691E+00 

19 

• i04E-02 

960.00 

•422229E+U2 

• 132142E+00 

20 

•iiOE-02 

909.4/ 

•  2  7024  /E  +  02 

.1223936+00 

2  1 

•ll6t-02 

864.00 

• 18 lb26E+u2 

•664643E-0i 

22 

•i22E-02 

822.86 

•  6  74702E  +  U1 

•  132660E+00 

23 

127E-02 

785.46 

•640139t+01 

•  283  73ot+00 

24 

.i33E-02 

731.30 

.103808E+U2 

.  1388  7  7t+00 

2b 

.  139E-02 

720.00 

•  1  1 228  7 E  +  02 

•  139  786E+00 

26 

•  14SE-02 

691.20 

.90  729  3E  +  01 

•2273396+00 

27 

.lSOt-02 

664.62 

• iu6933E+02 

•3167416+00 

.lb6E-U2 

040.00 

•  82  7  762E+01 

•3117636+00 

29 

.i62E-02 

61 7.14 

•  3  70629E+O1 

•3603396-01 

.16BE-02 

393.86 

•37131 7E+00 

.3028706-01 

31 

•  1  /‘♦t-02 

3  76.00 

•  89869  7t+U0 

•  33  72186-01 

TABLE  26 

FREQUENCY  RESPONSE  OF  THE  EARTH  AT  TBILISI  2  SEPTEMBER  1957 
(LOW-PASS  FILTERED). 
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FIGURE  38 

THE  FREQUENCY  RESPONSE  OF  THE  EARTH  VERSUS 
FREQUENCY  DETERMINED  FROM  EAST -WEST 
COMPONENT  OF  TELLURIC  FIELD  AT  TBILISI 
1  AND  2  SEPTEMBER  1957 
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FIGURE  39 

THE  FREQUENCY  RESPONSE  OF  THE  EARTH  VERSUS 
FREQUENCY  DETERMINED  FROM  NORTH-SOUTH 
COMPONENT  OF  TELLURIC  FIELD  AT  TBILISI 
I  AND  2  SEPTEMBER  1957 
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FIGURE  41 

THE  POWER  SPECTRA  OF  THE  NORTH- SOUTH 
MAGNETIC  AND  EAST-WEST  TELLURIC  FIELDS 
AT  TBILISI  2  SEPTEMBER  1957 
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TABLE  27 


COMPARISON  BETWEEN  FREE  EARTH  PERIODS  OBSERVED  WITH 
A  GRAVIMETER  AND  PEAKS  IN  THE  GEOMAGNETIC  POWER  SPECTRA 

Periods  are  for  the  spheroidal  mode  S. 

O  X/ 


Order, 

1 

Gravimeter 

Period 

Minutes 

E-W  Telluric  Period 

Minutes  ^  Difference 

N-S  Magnetic  Period 

Minutes  ’jo  Difference 

2 

54.98 

60.4 

-9.9 

60.4 

-9.9 

5 

35.87 

33.88 

+5.5 

38.4 

+7.05 

4 

25.85 

26.18 

-1.3 

26.18 

-1.3 

5 

19.83 

20.56 

-0.57 

20.56 

-0.57 

6 

16.07 

15.57 

+3.1 

16.00 

-0.4 

7 

15.42 

13.72 

-2.3 

14.05 

-4.7 

8 

11.78 

12.00 

-1.9 

12.26 

-4.1 

9 

10.57 

10.57 

+1.0 

10.86 

-2.7 

10 

9.685 

9.76 

-0.85 

no 

1  value 
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computed  power  spectra.  The  strong  components  (N-S  magnetic  and  E-W  telluric) 
on  2  September  were  chosen  for  comparison  since  these  data  are  the  most 
reliable  (see  Chapter  Vl). 

The  frequency  resolution  and  confidence  intervals  of  the  spectra 
presented  here  do  not  permit  positive  statements  about  this  interesting 
comparison.  It  is  only  suggested  that  physically  the  induced  currents  are 
possible  and  that  their  existence  may  be  evidenced  by  peaks  in  the  power 
spectrum  of  the  magnetotelluric  field. 
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